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PRAKATA 


Tugas penulis di PT. IPTN adalah mengembangkan teknologi Aerodinamika dan 
Dinamika Fluida pada umumnya, untuk menunjang kebutuhan perusahaan. Salah satu 
cabang penting dari teknologi yang harus dikembangkan adalah Dinamika Fluida 
Komputasional. Dalam melaksanakan tugas tersebut penulis dibantu oleh banyak staff, 
yang karena situasi dan kondisi yang ada direkrut dari sarjana-sarjana S1 yang baru lulus. 
Seorang staff pengembang Dinamika Fluida Komputasional (DFK) harus memiliki 
pengetahuan yang mendalam di bidang pembuatan program komputer, matematika, 
metoda numerik dan dinamika fluida. Pengetahuan ini, sampai tingkat yang memadai, pada 
umumnya belum dimiliki oleh seorang sarjana S-1 yang baru lulus. Namun demikian para 
staff baru tersebut tetap dituntut untuk melaksanakan tugasnya, baik sebagai pengguna 
maupun pembuat perangkat lunak DFK. Tentu saja ini menimbulkan masalah karena 
pengetahuan yang dimiliki belum cukup untuk menunjang tugas yang diemban. 

Disamping itu sejak tahun 1984 penulis telah berpengalaman menjadi pembimbing dari 
mahasiswa-mahasiswa yang ingin melakukan kerja praktek atau tugas akhirnya di PT. 
IPTN. Banyak di antara mereka akhirnya direkrut menjadi staff PT. IPTN, spesialis bidang 
aerodinamika atau dinamika fluida. Mahasiswa tersebut datang dari berbagai macam 
perguruan tinggi seperti ITB, ITS, ITI, UI, UGM, UNIBRAW, UNPAD dan bahkan dari 
Jerman. Latar belakang pendidikan mereka juga bervariasi seperti Matematika, Fisika, 
Teknik Mesin dan Teknik Penerbangan. Para mahasiswa dan staff baru tersebut pada 
umumnya mempunyai motivasi yang tinggi untuk belajar, tetapi sering mengalami 
kesulitan untuk belajar sendiri membaca buku-buku teks di bidang dinamika fluida 
komputasional. Masalah utamanya adalah karena tidak ada satu buku teks yang mencakup 
seluruh materi yang mereka butuhkan supaya segera dapat diterapkan dalam melaksanakan 
tugas. Informasi yang dibutuhkan tersebar di banyak buku teks dan jurnal-jurnal ilmiah, 
dan sering ditulis untuk pembaca yang berpengalaman sehingga sulit dimengerti oleh 
seorang pemula. Untuk membantu mereka, penulis terpaksa harus mengumpulkan, 
mengolah mempersiapkan bahan-bahan dalam bentuk tulisan tangan mengenai topik-topik 
yang dibutuhkan mereka. Diharapkan bahwa catatan-catatan tersebut dapat mempercepat 
dan memperlancar proses belajar para anak didik penulis. Catatan dan bahan-bahan yang 
dikumpulkan selama bertahun-tahun tersebut adalah bahan utama buku ini, setelah diolah 
dan diperbaiki lebih lanjut. Banyak bahan yang tadinya masih berwujud konsep kemudian 
akhirnya diterapkan oleh para anak didik menjadi program komputer yang dapat 
dimanfaatkan. Dengan demikian dapat pula dikatakan bahwa isi buku ini adalah rekaman 


dari sebagian hasil penelitian dan pengembangan yang dilakukan penulis dan dibantu oleh 
para anak didik tadi. Kepada semua mantan anak didik, yang banyak di antaranya sekarang 
menjadi kolega penulis, dengan ini penulis mengucapkan banyak terima kasih. Namun 
demikian segala kekurangan dan ketidak-sempurnaan yang ada dalam buku ini adalah 
tanggung jawab penulis sendiri. 

Harapan penulis adalah semoga buku ini bermanfaat untuk semua pembaca, yang di 
antaranya mungkin ada yang bercita-cita untuk menjadi seorang spesialis dinamika fluida. 
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BABI 
PENDAHULUAN 


Bagi seorang pemula masalah dinamika fluida atau pergerakan fluida, merupakan 
suatu masalah yang sangat pelik dan rumit karena fluida itu sangat mudah berubah bentuk 
dan pergerakannya sulit dibayangkan atau dijelaskan. 

Pergerakan solida (benda padat), seperti mobil, dapat dengan mudah dibayangkan. 
Mobil memiliki bentuk tertentu dan konsep pergerakannya, baik translasi maupun rotasi, 
dapat diterima akal tanpa kesulitan. Sebaliknya pergerakan fluida, misalnya aliran air 
sungai, itu sangat sulit dibayangkan apalagi dijabarkan. Semua benda memang terdiri dari 
partikel yaitu atom dan molekul. Tetapi partikel dasar tersebut begitu kecil dan berjumlah 
terlalu banyak, sehingga mustahil bagi kita- untuk mencoba menjelaskan aliran air sungai 
dengan menganalisa pergerakan masing-masing partikel air. Oleh karena itu yang dapat 
dilakukan adalah menganalisa pergerakan rerata statistis molekul-molekul yang berjumlah 
sangat besar itu. Inilah dasar pemikiran dari anggapan materi kontinyua, dimana fluida 
dianggap sebagai materi kontinyua, yaitu tanpa bentuk tetap seperti agar-agar dan bukan 
terdiri dari partikel-partikel diskrit seperti butiran pasir di pantai. 

Walaupun anggapan materi kontinyua itu memang memudahkan penjabaran 
pergerakan fluida, tetapi seperti telah dikatakan sebelumnya kita terbentur pada masalah 
konseptual yaitu bagaimana cara menjelaskan pergerakan benda yang bentuknya selalu 
berubah, tanpa batasan-batasan tertentu yang memberinya identitas tersendiri. Untuk 
mengatasi masalah tersebut ada 2 pendekatan yang dapat dilakukan yaitu pendekatan 
Lagrange dan pendekatan Euler. Dalam pendekatan Lagrange, kita membutuhkan konsep 
partikel fluida dengan identitas tertentu. Materi kontinyua atau agar-agar yang tanpa 
bentuk itu kita bagi menjadi sejumlah besar partikel-partikel fluida dengan identitas bentuk 
masing-masing dan saling berdempetan. Sebagai contoh, air dalam pipa dapat dibayangkan 
terdiri dari partikel-partikel air yang masing-masing berbentuk silinder kecil, misalnya 
berukuran panjang 1 mm dengan diameter 0.1 mm. Salah satu partikel kemudian dipilih 
untuk dipelajari pergerakannya. Untuk membedakan partikel tersebut dari partikel-partikel 
yang lainnya, kita bayangkan partikel yang dipilih itu diberi warna merah. Dalam 
pendekatan Lagrange, pergerakan partikel merah itulah yang diselidiki dan dianalisa, yaitu 
kecepatan gerak, percepatan, gaya-gaya yang beraksi dan perubahan bentuknya dianalisa 
sebagai fungsi waktu dan koordinat titik dimana dia berada. Untuk kasus aliran laminar, 
partikel merah tersebut akan mengalami pergerakan translasi dan rotasi, bentuk dan 
volumenya mengalami perubahan (deformasi), tetapi identitasnya tetap ada. Sebaliknya 


dalam aliran turbulen, identitas partikel merah itu pada akhirnya akan hilang. Dia akan 
terpotong-potong dan bercampur aduk dengan partikel-partikel lainnya, membentuk 
partikel baru. Kalau dalam aliran laminar identitas atau batasan partikel merah itu akan 
selalu tetap ada, sebaliknya dalam aliran turbulen suatu partikel fluida dapat "dilahirkan" 
atau terbentuk dan akhirnya "musnah", yaitu mempunyai "masa hidup" tertentu. Itulah 
sebabnya mengapa uraian tentang aliran turbulen menjadi jauh lebih rumit dari untuk kasus 
aliran laminar. 

Dari penjelasan di atas dapat dimengerti bahwa pendekatan Lagrange kita mengikuti 
pergerakan suatu partikel fluida tertentu dan mengamati sejarah pergerakannya sebagai 
fungsi koordinat ruang dan waktu. Perlu diperhatikan bahwa bentuk partikel fluida boleh 
dipilih sembarang, tidak terbatas pada bentuk silinderis seperti dalam contoh di atas. 
Ukuran besar partikel juga boleh diperkecil dengan catatan harus jauh lebih besar daripada 
ukuran molekular, sedemikian rupa sehingga dalam partikel fluida itu masih terdapat 
bermilyar-milyar molekul. Dari teori atom modern diketahui bahwa untuk gas dengan 
temperatur dan tekanan standar suatu volume sebesar 22.415 m?' mengandung sebanyak 
6.0248 x 10?" molekul (lihat misalnya Wehr and Richardo (1964)). Ini berarti bahwa suatu 
partikel fluida dengan volume 0.001 mm? , yaitu kubus dengan panjang sisi 0.1 mm, 
mengandung sekitar 2.7 x 10” atau 27 triliun molekul, suatu jumlah yang masih sangat 
besar. 

Dalam pendekatan Euler, kita tidak mengikuti pergerakan suatu partikel fluida 
tertentu. Sebaliknya kita memusatkan perhatian kita pada suatu volume kecil dengan 
bentuk tertentu, yang terpancang di titik tertentu dan disebut volume kontrol. Medan 
aliran dibagi menjadi sejumlah besar volume kontrol yang saling berdempetan. Salah satu 
volume kecil tersebut dipilih dan untuk memberinya identitas, kita bayangkan bahwa 
dinding volume kontrol ini diberi warna merah. Kecepatan dan gaya-gaya yang beraksi 
pada fluida yang berada dalam volume merah tersebut dianalisa sebagai fungsi koordinat 
ruang dan waktu. Apabila volume kontrol merah tersebut diperkecil, maka pada limitnya 
kita dapat berbicara tentang kecepatan dan gaya-gaya yang beraksi pada fluida di suatu 
titik sebagai fungsi waktu. Pergerakan fluida secara keseluruhan dapat dijabarkan dengan 
menghitung dan mengetahui kecepatan dan gaya-gaya yang beraksi pada fluida di semua 
titik dalam medan aliran sebagai fungsi waktu. Tentu saja titik-titik dalam medan solusi 
berjumlah tak-berhingga, sehingga dalam prakteknya kita hanya menganalisa pergerakan 
fluida disejumlah besar, tetapi berhingga, titik-titik tertentu dalam ruang dimana fluida 
berada, yang kita sebut sebagai medan solusi. Dari segi pengertian dan argumentasi fisis, 
sebenarnya pendekatan Lagrange itu lebih mudah diterima akal, karena lebih mirip dengan 
pergerakan solida seperti mobil. Tetapi secara matematis pendekatan lagrange menjadi 


sangat rumit dan jauh lebih sulit daripada pendekatan Euler. Inilah sebabnya mengapa 
pendekatan Euler lebih umum diterapkan dalam masalah pergerakan fluida. Pendekatan 
Lagrange hampir tidak pernah diterapkan kecuali untuk kasus-kasus tertentu dalam 
masalah aliran turbulen. Dengan bantuan konsep-konsep tersebut di atas maka upaya 
untuk menjelaskan pergerakan fluida dapat mulai dikembangkan. Penjelasan yang telah 
diberikan diharapkan dapat membantu para pemula yang ingin meniti kariernya sebagai 
pekerja atau peneliti di bidang dinamika fluida, untuk lebih menguasai permasalahan yang 
dihadapi. 

Masalah utama kedua yang menghadang para pemula di bidang dinamika fluida 
adalah kenyataan bahwa matematika yang terlibat dalam analisis dinamika fluida itu 
ternyata sangat sulit. Sebagai contoh, persamaan atur yang mengatur pergerakan fluida 
secara umum adalah persamaan Navier-Stokes. Sistem persamaan Navier-Stokes adalah 
suatu sistem persamaan diferensial parsial derajat 2, dalam empat (4) variabel bebas yang 
bersifat sangat tak-linier. Solusi persamaan ini sangat tergantung pada kondisi batas yang 
harus dipenuhi, sedangkan kondisi batas ini biasanya adalah benda 3-dimensional dengan 
geometri atau bentuk yang sangat rumit. Solusi analitis dari persamaan Navier-Stokes 
dapat diperoleh hanya untuk kasus geometri atau kondisi batas yang sangat sederhana saja. 
Untuk kondisi batas yang sedikit lebih rumit, solusi persamaan hanya dapat diperoleh 
secara numerik. Oleh karena itu para pemula di bidang dinamika fluida juga harus 
mempelajari matematika numerik supaya dia dapat melakukan analisis di bidang dinamika 
fluida. Kesulitan matematis tersebut di atas akan mengecilkan hati seorang pemula bila dari 
awalnya dia sudah disodori masalah dinamika fluida yang rumit. Untuk menghindari 
masalah tersebut dalam tulisan ini kita hanya akan membahas masalah dinamika fluida yang 
paling sederhana. Namun demikian masalah yang dibahas tetap harus bermanfaat secara 
praktis dan mampu memberikan landasan yang kokoh untuk mempelajari masalah-masalah 
selanjutnya yang lebih rumit. Tujuan dari buku ini adalah memperkenalkan beberapa 
masalah dinamika fluida, khususnya aerodinamika yaitu dimana fluidanya terbatas pada 
udara, kepada para pemula yang mungkin adalah sarjana fisika, matematika atau disiplin 
ilmu yang lain, dan tertarik untuk mempelajari aerodinamika teoritis secara mendalam, 
tetapi tidak memiliki latar belakang aerodinamika. Walaupun materi yang dibahas sangat 
dibatasi, tetapi diharapkan agar para pemula akan belajar banyak tentang konsep-konsep 
yang berlaku di bidang aerodinamika pesawat terbang. Para pemula juga diharapkan akan 
belajar banyak tentang konsep-konsep penyelesaian masalah persamaan diferensial baik 
secara analitis maupun komputasional. Pembaca akan diperkenalkan dengan konsep gaya 
angkat, gaya tahan, sirkulasi, vortisitas, garis arus, koefisien tekanan dan lain sebagainya 
yang sangat penting dalam aerodinamika pesawat terbang. Konsep-konsep tersebut juga 


diperlukan sebagai landasan untuk mempelajari masalah aerodinamika yang lebih maju dan 
jauh lebih rumit Teori variabel kompleks dan penerapannya dalam penyelesaian masalah 
aliran potensial tak-mampat juga dibahas. Pemetaan konformal yang menghasilkan bentuk 
aerofoil dibahas dengan tuntas dan digunakan untuk mempelajari parameter geometri 
aerofoil, yaitu ketebalan, kelengkungan dan kelancipan ekor, dan pengaruhnya pada sifat- 
sifat aerodinamika aerofoil. Penyelesaian masalah aliran potensial dengan metoda panel 
dan metoda taksiran beda berhingga juga dibahas sebagai pengenalan awal dengan masalah 
aerodinamika komputasional pada umumnya. 


BAB II 
PEMETAAN KONFORMAL DENGAN 
METODA DERET LAURENT 


Dari teori variabel kompleks kita tahu bahwa solusi dari persamaan Laplace 2- 
dimensional adalah semua fungsi analitik kompleks. Model matematis dari masalah aliran 
tunak, potensial, tak-mampat dan 2-dimensional adalah persamaan Laplace, jadi masalah 
aliran tersebut secara teoritis dapat diselesaikan dengan bantuan teori variabel kompleks. 
Sebagai contoh, aliran potensial di sekeliling bentuk lingkaran dapat dianalisa secara 
tuntas karena rumus analitis untuk fungsi potensial aliran ini dapat dicari dengan mudah. 
Sebaliknya fungsi potensial kompleks untuk aliran di sekeliling aerofoil sembarang itu 
tidak dapat dicari secara langsung. Namun demikian diketahui bahwa bentuk aerofoil itu 
dapat dipetakan menjadi bentuk lingkaran oleh suatu rumus pemeta berbentuk fungsi 
analitik. Disamping itu diketahui pula bahwa semua fungsi dari fungsi analitik pasti juga 
bersifat analitik dan oleh karena itu juga merupakan solusi persamaan Laplace. Apabila 
rumus pemeta lingkaran menjadi aerofoil yang diinginkan itu diketahui, maka rumus 
fungsi potensial kompleks untuk aliran di sekeliling aerofoil tersebut dapat dicari secara 
tidak langsung dan masalah dapat diselesaikan. Ini berarti bahwa masalah analisis aliran 
potensial di sekeliling bentuk aerofoil sembarang dapat disederhanakan menjadi masalah 
mencari rumus pemeta bentuk lingkaran menjadi bentuk aerofoil tersebut. Dalam hal ini 
kita akan membahas masalah pemetaan tersebut. 

Dari teori variabel kompleks diketahui bahwa di kawasan cincin ("annulus") yang 
dibatasi oleh 2 lingkaran konsentris di bidang kompleks C , yaitu di kawasan 


Ratam s Ic 5 Gol & Rar 
suatu fungsi analitik kompleks Z(C) sembarang itu selalu dapat diwakili oleh suatu deret 


Laurent berikut 


Ha ap " e.1) 


m-—1 


dimana untuk m-—0,1,2,... 


CA Ya 
"Ima CG)" 
Lg anjar 


"ami! (CG) 
Bukti untuk pernyataan diatas dapat dilihat di buku teks tentang variabel kompleks dan 


(2.2) 


juga telah dibahas oleh Wilianto (1992), jadi tidak perlu dibahas disini. Tanpa membatasi 
manfaat hasil pembahasan dan untuk mempermudah penulisan, kita anggap bahwa C, 
berada di titik asal sistem koordinat yaitu G, - (0,0). 

Pola aliran angin seragam dikejauhan tak-berhingga dari lingkaran atau aerofoil 
itu tidak terganggu oleh keberadaan lingkaran atau aerofoil. Ini berarti bahwa dikejauhan 
tak-berhingga z — G atau 

Tini z5 (2.3) 


Persyaratan diatas akan terpenuhi apabila 


B,,—0 untuk m-2, 3, 4,... | C1) 


B, 51 untuk m-1 
Penerapan persyaratan diatas pada rumus (2.1) dan (2.2) akan memberikan rumus 
berikut sebagai rumus pemeta yang diinginkan. 


Z6) 6 Ya (2.5) 
C, - 8: page” ag untuk m- 0, 1, 3... (2.6) 
21i 


Karena deret di ruas kanan rumus (2.5) itu konvergen, maka suku-suku wa 


diperkirakan semakin mendekati nol untuk harga m yang semakin besar. Ini berarti 
bahwa deret tersebut boleh dipancung setelah suku ke M dengan galat yang kecil tak 
berarti, apabila M bernilai cukup besar. Oleh karena itu untuk tujuan hitungan praktis, 
deret (2.5) dapat dipancung menjadi 


A0 La Fa (27 


Masalah yang dihadapi sekarang adalah bagaimana cara menghitung nilai koefisien- 
koefisien kompleks C,,. Misalkan bentuk aerofoil yang ingin dianalisa itu diberikan, yaitu 
sebagai kurva ZA di bidang-z. Spesifikasi bentuk aerofoil diberikan sebagai sejumlahan 
koordinat titik-titik, misalkan ZAA)) untuk j — 1, 2, 3,...,JTA, yang berada pada 
permukaan aerofoil tersebut. Untuk pembahasan selanjutnya ditentukan bahwa aerofoil 
harus berekor runcing, yaitu titik tepi buritan atas dan bawah harus berhimpit. Kita 
tentukan juga bahwa titik buritan atas adalah titik ZAA(1) sedangkan ZAA(JTA) adalah 
titik tepi buritan bawah dan nilai indeks j meningkat arah kebalikan putaran jarum jam. 
Aerofoil ZA dipetakan oleh rumus pemeta (2.7) menjadi lingkaran berjari-jari R di 
bidang-C, yaitu lingkaran GL dimana (cu R. Secara lebih khusus titik ZAA() 
dipetakan menjadi titik (LA (j) yang berada pada lingkaran CL. 

Walaupun koordinat ZAA(j) diberikan tetapi koordinat titik petanya, yaitu 
CLA(j), itu tidak diketahui jadi rumus (2.6) tidak dapat diterapkan secara langsung 


untuk menghitung nilai koefiasien-koefisien deret Laurent terpancung (2.7). Koordinat 
titik lingkaran CLA (j) adalah (R,0A(j)) dan baik R maupun distribusi sudut OA) itu 


tidak diketahui. Untuk sedikit lebih mempermudah masalah, kita rumuskan suatu rumus 


transformasi baru sebagai berikut 


SA 23) 


Rumus (2.8) dapat memetakan lingkaran feri 5 R menjadi lingkaran satuan, berjari-jari 
1, di bidang-£” yaitu lingkaran GS dimana (CS| - 1. 


Rumus deret Laurent terpancung (2.7) dan rumus untuk koefisiennya (2.6) 
sekarang dapat ditulis ulang menjadi 





Ac) -Rc 3 m (2.9) 
“s(c) 
“Em Ne) ae 
aa Aa dc)" ac (2.10) 


Peta dari titik ZAA) adalah titik CSA(j) dibidang &“ dengan koordinat (1, OA(j)). 


Pada umumnya @A(1) tidak bernilai nol dan ditentukan bahwa @A(1)— A yang belum 


diketahui nilainya. Untuk memudahkan hitungan, kita menginginkan sudut yang pertama 


tersebut bernilai nol dan untuk tujuan tersebut ditentukan bahwa 
GC" —6' exp(-ih) (2.11) 
Peta dari titik CSA(j) dibidang €” adalah titik CSP(j) dengan koordinat (1, 6P/))) 


dimana 

OP(j) -@A(j)-A (2.12) 
dan khususnya 

@P(1)-—0 (2.13) 


Rumus pemetaan (2.9) dan (2.10) sekarang dapat ditulis ulang menjadi 





ic”) — R exp(ih)C” 4 » Cm 2.14) 
“(c”) 
Ma Unt 201 pen oo 
Ch - exp(im2) P —a pal ytc ) dc (2.15) 


Untuk memudahkan penulisan tanda bintang akan dilepas, jadi rumus pemetaan yang 
mampu memetakan aerofoil ZA di bidang-z menjadi lingkaran satuan CS di bidang-£ 
dengan catatan bahwa titik ekor aerofoil ZAA(l) dipetakan menjadi titik lingkaran 
dengan koordinat (&, ep(1)) dimana @P/(1) — 0, sekarang ditulis ulang dari persamaan 
(2.14) dan (2.15) sebagai berikut 


XC) — R exp(iA)t 4 re (2.16) 
1 m4 
San hao)er tag (2.17) 


Titik ZAA(j) pada permukaan aerofoil ZA di bidang-z dipetakan oleh rumus-rumus 
(2.16) dan (2.17) menjadi titik-titik (1, OP(j)) pada lingkaran |GSP| -1 di bidang-C 
dengan pengertian bahwa @P/1) selalu bernilai nol. Karena kurva ZA dipetakan menjadi 


lingkaran (SP maka persamaan (2.17) dapat disederhanakan menjadi 


Cm -—— $ZA . CSP" dCSP 
2ni 


Karena GSP adalah lingkaran satuan pada bidang-C maka titik-titik pada lingkaran 


tersebut mempunyai koordinat (1, 0) dan rumus diatas dapat ditulis menjadi 


Oa Pa pzA . explilm-1)0) alexp (i0)) 
atau 
C,- 5-pza . exp(imO) do (2.18) 


Data asli bentuk aerofoil memang diberikan dalam bentuk sejumlah terbatas koordinat 
titik-titik ZAA untuk 1S jX JTA. Suatu rumus interpolasi dapat diturunkan untuk 
mewakili bentuk aerofoil dengan menggunakan data titik-titik ZAA(I) ke ZAA) 
tersebut. Panjang kurva dari titik ZAA(1) ke ZAAG) sepanjang permukaan aerofoil dapat 
dihitung dan diberi simbol SAA(j), misalnya. Komponen real dan imajiner ZAA(j) adalah 
XAA(j) dan YAA). Suatu fungsi interpolasi x(s) yang melewati himpunan pasangan 
data (SAAG), XAAG), j — 1, JTA) dapat dihitung. Hal yang sama dapat dilakukan untuk 
mencari rumus interpolasi y(s) yang melewati semua pasangan (SAA), XAAM)). 
Apabila harga s diberikan maka koordinat titik pada permukaan aerofoil dapat dihitung 
dengan rumus-rumus interpolasi x(s) dan y(s) tadi. Dengan cara tersebut semua titik pada 
permukaan aerofoil dapat dihitung koordinatnya sebagai hasil prosedur interpolasi dan 
diberi simbol ZAI, dengan koordinat (XAI, YAI) untuk panjang kurva permukaan 
aerofoil dapat dihitung koordinatnya sebagai hasil proses interpolasi, sekarang masalah 
bagaimana cara menghitung koefisien C,, dengan rumus (2.18) dapat dijelaskan sebagai 
berikut 

Lingkaran satuan CSP ditentukan diwakili sejumlah, JTI, titik-titik yang 


terdistribusikan secara seragam. Koordinat titik-titik tersebut adalah (&, ep(i)) dimana 


@P(1)-0 





@P(i--1)-i. untuk i — 1,2, ...,JTI-1 


JTI-1 
Peta dari titik GSP(i) adalah titik ZAI(i) pada permukaan aerofoil ZA di bidang-z. 
Koordinat titik ZAI(i) pada umumnya tidak diketahui kecuali untuk i—1 dan i-JTI karena 
kedua titik tersebut selalu berhimpit dengan titik ekor (tepi buritan) aerofoil. Namun 
demikian diketahui bahwa ZAI berada pada permukaan aerofoil ZA. 


Pada awalnya distribusi panjang kurva SAI) dari titik ZANI) ke ZANIi) itu 
ditebak, misalnya dengan menganggap bahwa titik-titik ZAI(i) itu terdistribusi secara 
seragam sepanjang permukaan aerofoil, yaitu 

SAI) —0 

SAK 1) - SANi) -ASAI untuk i-1,2,..., JTI-2 
Nilai SAI(JTI) itu sudah diketahui karena ini adalah panjang seluruh permukaan aerofoil 
yang dapat dihitung dari data masukan, dan SAI(JTI) - SAA(JTA) -— SSA, misalnya. 
Nilai ASAI dapat dihitung sebagai berikut 


- SSA 
ASAT A. IT» 


Rumus interpolasi x(s) dan y(s) kemudian diterapkan untuk menghitung koordinat titik- 
titik ZAIG), yaitu (KAN), YANi)) untuk semua nilai i. 


Karena koordinat ZAI dan petanya, yaitu (4, or(i)), sekarang diketahui maka 


rumus (2.18) dapat diterapkan untuk menghitung semua nilai koefisien C,,. Tentu saja 
nilai C,, yang diperoleh adalah nilai tebakan karena ZAI yang digunakan dalam hitungan 
adalah nilai koordinat tebakan. Untuk mempertegas kenyataan tersebut nilai C,, tebakan 
tadi diberi simbol CT,,. Pemasukan nilai-nilai itu ke dalam persamaan (2.16) akan 


memberikan koordinat titik-titik pada arofoil tebakan ZTI, yaitu 


ZTI —R exp(i(A #OP))- 5 CT,, exp(-imOP) (2.19) 


m-0 
Koordinat ZTI(i) adalah tebakan dari koordinat ZAI(i). Apabila setiap ZTI(i) bernilai 
sama dengan ZAI(i) maka tebakan yang dilakukan sudah tepat dan CT, adalah nilai dari 
C,, yang diinginkan. Bila ZTI(i) belum sama dengan ZAN(i) untuk semua i maka 
distribusi ZAI(i) disepanjang permukaan aerofoil ZA harus ditebak ulang atau harus 
diperbaiki sebagai berikut. 

Walaupun aerofoil tebakan ZTI berbeda bentuk dari ZA yang diinginkan, tetapi 
distribusi panjang STI(1) dari ZTN1) ke ZTN(i) adalah taksiran yang lebih baik dari 
distribusi SAI(i) yang sebenarnya bila dibandingkan tebakan distribusi SAN) yang 
digunakan sebelumnya. Oleh karena itu tebakan untuk SAI(i) kemudian diperbaiki 
sebagai berikut 


SSA 


(2.20) 
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Tindakan selanjutnya adalah mengulang seluruh prosedur yang telah dijelaskan 
sebelumnya sampai akhirnya proses iterasi (tebak ulang) menjadi konvergen, yaitu semua 
ZTX(i) menjadi sama dengan ZAN(). 

Dalam prosedur di atas kita menganggap bahwa nilai untuk R dan 4 dalam rumus 
(2.19) itu telah diketahui, sedangkan dalam kenyataannya nilai-nilai tersebut belum 
diketahui. Nilai R dan 4 untuk setiap iterasi dapat dihitung sebagai berikut. 

Dalam bab 7 modul I bahan kursus ini telah ditunjukkan bahwa titik ekor aerofoil 
adalah titik kritis yang harus memenuhi persyaratan berikut 





dz 

—20 2.21 

.. 2.21) 
Turunan dari rumus pemeta (2.16) adalah 

dz Cc 

— SR exp(i2)-) m — 2.22 

ae Bapa) Ym Ta (2.22) 


Pemasukan (2.22) ke (2.21) memberikan hasil berikut untuk G yang mewakili peta dari 


titik ekor aerofoil 
M 
Rexp(i4) — 9 mC,, (2.23) 
ml 
karena GSP(i) — (1, 0). 


Persamaan (2.23) menunjukkan bahwa R adalah modulo dan 4 adalah argumen 


M 
dari konstanta kompleks Y'm C,, . Karena tebakan nilai C,,, yaitu CT,,, itu diketahui 


m-l 


maka rumus (2.23) dapat digunakan untuk menghitung R dan 4 yaitu 


R- mod 3, cx,) 


m-1 


aa Xon,| 


m-i 


(2.24) 


Dalam metoda yang dijelaskan di atas, aerofoil ZA di bidang-z secara langsung 
ditransformasikan menjadi lingkaran (csP| 51 di bidang-G. Disamping metoda langsung 


tersebut juga ada metoda tak-langsung yang dapat diterapkan. Pertama-tama aerofoil 
ditransformasikan menjadi bentuk hampir lingkaran dengan menggunakan rumus pemeta 
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yang sederhana seperti rumus Joukowski atau Karman-Trefftz. Setelah itu bentuk hampir 
lingkaran yang diperoleh dipetakan menjadi bentuk lingkaran dengan menggunakan 
rumus deret Laurent. 

Rumus pemeta Karman-Trefftz telah dibahas secara rinci dibab 7 modul I bahan 
kursus ini, jadi hanya perlu dibahas secara ringkas di bawah ini. Lingkaran 


|Z2L| -1 di bidang-z' dapat dipetakan menjadi aerofoil Karman-Trefftz adalah 


(z ta)" tra) 
(z ta)" Ar -a) 


z(z ) —ka 





(2.25) 


Sedangkan rumus balik pemeta aerofoil Karman-Trefftz menjadi lingkaran adalah 





(2 kat 4 (z—ka)/ 
Ag 
(z-- ka) — (z— ka) 


Rumus Joukowski dapat diperoleh dari rumus Karman-Trefftz di atas dengan 


(2.26) 


menggunakan nilai k-2, yaitu 


ZTN) — R exp(ih). exp(ieP/(i)) -- Ser, exp(-imOP/i)) 


m-0 


2 


zz) 5z 4 2 (2.27) 
z 





dan rumus baliknya adalah 


2) ata) (2.28) 


Apabila rumus balik Karman-Trefftz (2.26) diterapkan pada suatu bentuk aerofoil 
sembarang yang mirip dengan bentuk aerofoil Karman-Trefftz, maka bentuk hasil 
pemetaan tentu saja bukan lingkaran tetapi hampir lingkaran. 

Untuk kasus dimana ZB dan ZD adalah koordinat titik dari titik ekor dan hidung 
aerofoil yang berada pada sumbu-x, rumus (2.26) dapat disederhanakan menjadi 





ph ala-zB) 
-e-n)6 ZD/« #(z—ZB) (2.29) 


2k J (z-2D)4—(z—2B)« 


Parameter k adalah parameter kelancipan ekor aerofoil, yaitu 
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K2 (2.30) 
T 


dimana t adalah sudut (dalam satuan radian) yang diapit oleh tangen atas dan tangen 
bawah pada titik ekor aerofoil yang ingin dipetakan. 

Hasil pemetaan adalah kurva hampir lingkaran ZHL. Titik pusat luas ZHL adalah 
ZHP dengan koordinat yang dapat dihitung sebagai berikut 


za Y rana) (2.31) 


dimana ZHL(n) adalah koordinat titik-titik hampir lingkaran. Sebelum dipetakan menjadi 
lingkaran |CSP|—1 di bidang-C, bentuk hampir lingkaran ZHL harus digeser dulu 


sehingga titik pusat luasnya menjadi berhimpit dengan titik asal sistem koordinat. Untuk 
tujuan pembahasan selanjutnya kita menganggap bahwa titik pusat ZHL selalu berada di 
titik asal sistem koordinat. Disamping itu untuk memudahkan penulisan simbol bintang 
akan dilepas, yaitu hampir lingkaran ZHL berada di bidang-z dan dipetakan menjadi 
lingkaran CSP di bidang-C. 

Rumus pemeta yang digunakan adalah rumus deret Laurent terpancung (2.16) 
dan (2.17), sama dengan untuk kasus pemetaan langsung bentuk aerofoil menjadi 
lingkaran. Prosedur hitungan nilai-nilai koefisien C,, pada dasarnya juga sama seperti 
telah dijelaskan sebelumnya. Rumus-rumus interpolasi x(s) dan y(s) yang mewakili 
bentuk hampir lingkaran dapat dicari. Prosedur iterasi untuk menghitung C,, itu sama 
dengan untuk pemetaan langsung bentuk aerofoil, dengan satu perbedaan mendasar yaitu 
bahwa hampir lingkaran ZHL tidak memiliki titik kritis. Ini berarti bahwa kita tidak dapat 
menerapkan rumus (2.24) untuk menghitung R dan A. Nilai R dan 4 untuk kasus ini 
harus dihitung dengan cara lain sebagai berikut. Setelah nilai C,, tebakan, yaitu CT.,, 
dihitung maka kita dapat menggunakan rumus untuk kurva tebakan (2.19) sebagai 
berikut 


ZTI(i) —R exp(ih). exp(iOP(i)) Ser, expl-imeP(i)) 


m-0 


atau 


R exp(ih) — (ea 2 » CT, eol-imor(G) |exo-ior6) 


m-—0 
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Apabila proses iterasi sudah konvergen maka pada akhirnya ZTI(i) - ZHI(i) dimana 
ZHI(i) adalah koordinat titik pada kurva ZHL sebagai tebakan titik peta dari titik 


lingkaran CSP(i) dan dihitung dengan rumus interpolasi seperti telah dijelaskan 
sebelumnya. Walaupun proses iterasi belum konvergen nilai ZTI(i) tidak akan berbeda 
jauh dari nilai ZHIi). Oleh karena itu rumus untuk R. exp(iA) di atas dapat ditaksir 
sebagai berikut 


R exp(ih) « (am - » CT, ep-imor()) lesn-ior6) (2.32) 


m-0 
Untuk memperkecil kesalahan yang terlibat dalam taksiran diatas, kita dapat mengambil 
nilai rata-rata berikut 





MN. £ (fana cr, sotimor()) Lesson) 


JTI-1 45 m0 
(2.33) 
Dari rumus (2.18) diketahui bahwa 


1 
C, - —?ZHI.dOP 
02x ? 
Integral di atas dapat dihitung secara numerik dengan metoda trapesium sebagai berikut 
1 JTI-1 
-— Y ZHIi).A0P 
Ti 


dimana 


- Mx 
AP ng 


jadi 





oJTI- 1 4 Yam 


Pemasukan nilai C, di atas ke rumus (2.33) akhirnya memberikan taksiran berikut 





Rexp/iA) ae Iran .exp(- fm or) (2.34) 


m-1 
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Rumus (2.34) dapat digunakan untuk menghitung nilai-nilai taksiran R dan 2 dalam 
kasus metoda tak-langsung atau dalam kasus pemetaan hampir lingkaran menjadi 
lingkaran. 

Rumus pemeta deret Laurent yang diperoleh adalah rumus (2.16) yaitu 


AO) - Rebana Ya Sa @19 


Rumus di atas dapat ditulis ulang menjadi bentuk yang diberikan oleh rumus (2.7) yaitu 


Mc 
Ao) ti (2.35) 
“3 (c) 
apabila ditentukan bahwa 
C"—R exp(iA)t (2.36) 
C,, -R" exp(imA)C,, (2.37) 


Perlu diperhatikan bahwa rumus (2.16) akan memetakan aerofoil ZA, yang 
diberikan sebagai data, menjadi lingkaran satuan dibidang-C dan titik ekor aerofoil 
dipetakan menjadi titik (1, 0) pada lingkaran tersebut. 

Disisi lain, rumus (2.35) akan memetakan aerofoil ZA menjadi lingkaran dengan 
jari-jari R dan titik ekor aerofoil dipetakan menjadi titik (R,A) pada lingkaran tersebut. 
Dalam prakteknya kita menginginkan proses pemetaan yang kedua, yaitu yang diberikan 
oleh rumus pemeta (2.35). Oleh karena itu setelah prosedur iterasi untuk menebak harga 
CT, itu konvergen, maka harga R dan A terakhir yang diperoleh harus disimpan dan 
rumus (2.37) digunakan untuk menghitung nilai C', yang dibutuhkan dalam rumus 
(2.35). Rumus pemeta (2.16) sesungguhnya hanya merupakan alat bantu dalam upaya 
menghitung nilai-nilai R, A dan koefisien-koefisien deret Laurent C,,. Dalam 
pembahasan selanjutnya kita akan menganggap bahwa rumus pemeta deret Laurent 
adalah rumus yang diberikan dalam bentuk (2.35) dan bukan dalam bentuk rumus (2.16), 
tetapi dengan catatan bahwa nilai R dan A telah dihitung dan disimpan. 
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BAB III 
METODA-METODA NUMERIK DALAM PENERAPAN 
PEMETAAN KONFORMAL DENGAN DERET LAURENT 


Konsep konsep yang sudah dijelaskan dalam bab sebelumnya dapat diterapkan 
untuk membuat program komputer untuk menganalisa karakteristik aerodinamik bentuk 
aerofoil sembarang. Prosedur hitungan yang perlu dilakukan dapat dijabarkan sebagai 
berikut. 

1. Geometri aerofoil ZA yang ingin dianalisa diberikan sebagai data masukan, yaitu 
sebagai sejumlah JTA koordinat titik-titik ZAAG) —- (XAAG) ,YAAG)) untuk 
1S jXJTA. Titik ZAA(I) berhimpit dengan titik ZAA(ITA) dan berada ditepi 
buritan aerofoil, sedangkan nilai indeks j meningkat kebalikan arah putaran jarum 
jam. 

2. Koordinat "titik pusat" aerofoil dihitung yaitu 


1 JTI-1 


ZAP — ZAA(j 
JTI-1 z @ 





Aerofoil kemudian digeser sehingga titik pusatnya menjadi berhimpit dengan titik 
asal bidang-z 
ZAA(j) - ZAA(j) -ZAP untuk 1£ j « JTI 

3. Distribusi jarak SAA(j) sepanjang permukaan aerofoil dari ZAA(1) ke ZAA) 
dihitung dan dinormalisasikan. 


SAAM) 
SAA(JTA) 
4. Rumus interpolasi x(s) melewati himpunan titik-titik data (SAAG), XAAG)) untuk 

1x js JTA ditentukan. 
5. Rumus interpolasi y(s) melewati himpunan titik-titik data (SAAG), YAAG)) untuk 
1X j«x JTA ditentukan. 


SAA() — untuk 1x js JTA 


6. Koordinat titik-titik CSP(i), yaitu (1,0P(i)) untuk 1x iX JTI, ditentukan sebagai 
berikut 
@P(1)-—0 


untuk 1Xis JTI-1 





Opfi-x1) 5ia — 7 
JTi 


11. 


13. 


Nilai awal panjang kurva SAIXi) sepanjang permukaan aerofoil ZA dari titik ZAI(1) 
ke ZAI(i) ditebak. ZAN1) berada ditepi buritan aerofoil sedangkan ZAI(i) adalah 
tebakan dari titik peta pada aerofoil ZA dari titik lingkaran CSP(i) 

SAX(i) —0 

i 
JTI— 

Koordinat titik ZAI(i), yaitu (KAKI) YANG), pada aerofoil ZA dihitung dengan 
menggunakan rumus interpolasi berdasarkan tebakan distribusi SAN). 

ZAI(1) — ZAA (1) 

ZANJTI) — ZAA (JTA) — ZAA (1) 
Interpolasi untuk ZAI(i) - (KAK), YAIG)) 





SAK 4-1) — 7 untuk Si JTI-1 


XAKi) - x(s - SAKi)) 
YAH) — y(s — SAK)) 


Tebakan nilai koefisien deret Laurent dihitung 


| untuk 2 Sis JTI-1 


CT(m) - 55 $ZAL.cep(imOP)dOP untuk 0s m «XM 


Tebakan nilai R dan 4 dihitung 


R.exp(iA) — Sim.CT(m) 


Koordinat titik ZTI pada aerofoil tebakan ZT, peta lingkaran GSP, yang diperoleh 
dari pemetaan dengan menggunakan rumus deret Laurent tebakan dihitung sebagai 
berikut. 


M 
ZTKi) - R.exp(i(a-OP())) « 9 CT(m).exp(-imeP/i)) 
m—0 
untuk1 Si JTI-1 
Galat koordinat hasil tebakan dihitung. 
egalat(i) - mod(ZTKi) -ZANi)) untuk 1 Si JTI-1 
Galat kemudian disortir untuk menentukan nilai galat terbesar yaitu galat max. 
Ketepatan hasil tebakan diuji. Apakah galatmak « & ?. Kuantitas & adalah 
bilangan positif kecil yang ditentukan nilainya, misalnya £ — 10“, sebagai ukuran 
ketepatan hasil tebakan. 
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14. 


Apabila jawaban dari pertanyaan di atas adalah tidak maka tebakan harus diperbaiki 
dengan memperbaiki distribusi panjang kurva SAK(). Distribusi panjang kurva ST(i) 
dari ZTX1) ke ZTI(i) sepanjang permukaan aerofoil ZT dihitung dan digunakan 
sebagai distribusi SAI(i) yang baru, yaitu 


ST) 
ST(JTI) 5 
selanjutnya seluruh prosedur diulang lagi atau diiterasikan mulai dari prosedur pada 
butir 8 
Apabila jawaban pada butir 13 adalah ya maka proses iterasi sudah konvergen dan 


SAK) — untuk 1ISis JTI 


nilai-nilai koefisien deret Laurent, R dan A yang benar sudah diperoleh. Koefisien 
deret Laurent CT(m) untuk rumus (2.16) sekarang diubah menjadi koefisien C(m) 
untuk rumus deret Laurent bentuk (2.35) yaitu 


C(m) - R".exp(imA).CT(m) untuk 0x ms M 


rumus pemeta deret Laurent adalah 


v(cLv) 
V(ZAV)-— ye ( av) 


Lingkaran lcl —R dipetakan oleh rumus di atas menjadi aerofoil ZA dengan catatan 
bahwa titik (R,A) dipetakan menjadi titik tepi buritan (ekor) aerofoil. 

Nilai sudut serang d dibaca sebagai data masukan. Koordinat titik stagnasi belakang 
lingkaran adalah (R,A) atau (R,-al) karena @ bernilai negatif. Kekuatan vorteks 
adalah y 41 sin|P| dimana |P| -a JA dengan anggapan bahwa @ selalu bernilai 


positif. Fungsi potensial kompleks untuk aliran ini adalah 








Oo) Cop ia) AR aptia) AR ac 
Cc 2x 
atau 
R?.explia) 





DG) - G.crp(—ia)-t 42R sinBI-In 


dimana € —(r,@) dan r 2 R. 


Bila dirumuskan bahwa 
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ear 


maka rumus di atas dapat ditulis menjadi 





R?.exp(ia) Ha 2R sin|B| 
E: Ng 


ea Sa MEN 20) 


O'(C) — exp(-ia) 





Rumus untuk cw) dapat ditulis dalam variabel r' dan 0 , sebagai berikut 





@(r”,0)-r”.exp(i(0 a))t--expl i(@ @))--2.sinBI(—6 4 imi”) 


Komponen-komponen riil dan imajiner dari fungsi potensial komples di atas adalah 





ae.0)-(r eh Joo x)-2.singI.0 


W(r",0) (F 1 Hino —)42.singl-In 


Turunan dari fungsi potensial kompleks di atas adalah 


@' (6) — exp(—ia) R?.exp(ia) |, 2Rsin|B| 








& € 
La SA 


Kurva isopotensial dapat dihitung koordinatnya dengan menyelesaikan persamaan 


6(r,@)-K, sedangkan garis arus dari persamaan w(r,0)-C dimana K dan C 


adalah konstanta yang diketahui. Kecepatan di titik (r,0) adalah 


lo'(r,0) - lo'(r',o), yang dapat dihitung sebagai berikut 
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V? -A?4B? 














Nagan cos(20 -“) sn 2sin|B sin 0 
ba (re)? r 

Ba sin(2 -a) da 2sin|P| cos@ 
Fa (r3) r 


Untuk menghitung koordinat titik-titik pada garis arus aliran di sekeliling aerofoil, 
prosedur yang harus dilakukan adalah 


16.1 Hitung koordinat titik-titik pada garis arus di sekeliling lingkaran el: 
dengan menyelesaikan persamaan w(r”,0) —K. Misalkan titik CGA" dengan 


koordinat (rca",0Ga) adalah titik pada garis arus tersebut. 
16.2 Garis arus di sekeliling garis arus ZA kemudian di hitung koordinatnya dengan 


memetakan titik-titik dengan koordinat (R.rGA",0GA) menggunakan rumus 


pemeta Z(C) —4 3 ae dimana € — (R.rGA',0GA). 


Koordinat kurva isopotensial di sekeliling bentuk aerofoil dapat diperoleh dengan 

cara yang sama. Untuk menghitung distribusi kecepatan aliran di sepanjang 

permukaan aerofoil ZA, prosedur yang harus dilakukan dapat dijabarkan sebagai 

berikut. 

16.3 Misalkan titik ZAV adalah titik pada permukaan aerofoil ZA dimana nilai 
kecepatan aliran ingin dihitung. Koordinat dari titik GLV yang merupakan 
peta dari ZAV dapat dihitung dari rumus transformasi deret Laurent z(C) 


yang koefisiennya telah dihitung. Misalkan koordinat CLV adalah 
(RLV,0LV) maka kecepatan di titik tersebut dapat dihitung dari rumus 


untuk lingkaran satuan yaitu V(CLV) — loe) 





dimana €' — (F8 0) 


16.4 Kecepatan aliran di titik ZAV kemudian dapat dihitung sebagai berikut 
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Va) 
V(ZAV) - Ye ( av) 


dimana 


z'(C) adalah turunan dari rumus pemeta, yaitu 


MO Aa Pm an 


Seperti dapat dilihat, bagian paling sulit dari masalah aliran di sekeliling aerofoil ZA 
adalah menghitung koefisien deret Laurent yang mampu memetakan aerofoil ZA 
menjadi bentuk lingkaran. Sekarang marilah kita pelajari secara lebih mendalam 
beberapa metoda numerik yang dibutuhkan dalam prosedur yang dijelaskan di atas. 
Untuk menghitung koordinat titik-titik pada garis arus di sekeliling lingkaran, kita 
harus menyelesaikan sebuah persamaan berbentuk sebagai berikut 


(—1Jsinto -a)- 2sinflinr — w (8.1) 

r 

dimana «,B dan w adalah kontanta-konstanta yang diketahui nilainya. Pola aliran 
yang diberikan oleh persamaan (3.1) pada dasarnya tidak dipengaruhi oleh sudut « 
karena @ disini hanya berfungsi sebagai fungsi pemutar. Apabila f(X max) «0 
adalah solusi dari persamaan (3.1) dimana &« -0, maka untuk & #0 solusinya 


adalah (RS,0S #a). Dengan pengertian tersebut persamaan (3.1) kita sederhanakan 





menjadi 
1. 
(e-#)sinosnnr-v (3.2) 
r 
dimana 
B-2sin|P| (konstanta) (8.3) 
Dimana nilai r kita tentukan, misalnya r —- RS, maka nilai 0 dapat dihitung dengan 
mudah yaitu 
— B.InRS).RS 
@ — sin" (w - (8.4) 
RS -1 
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dengan mengingat bahwa r » 1. Sesungguhnya nilai r terkecil yang diperbolehkan itu 


tergantung pada nilai w. Nilai r minimal adalah nilai r untuk 0-5 atau 


(w-B.Inr)r—r?—1. Apabila nilai r yang dipilih ternyata bernilai lebih kecil 
daripada nilai r minimal maka persamaan (3.4) tidak memiliki solusi. 
Untuk menentukan nilai r minimal tersebut, misalnya RMIN, kita harus mencari 


solusi dari persamaan f(r) — 0 dimana 

f(x) —r' —(w—B.Inr)r—1 (8.5) 
suatu cara yang lebih baik untuk menghitung kordinat titik-titik pada garis arus 
adalah dengan mengubah bentuk persamaan polar (3.2) menjadi persamaan dalam 
sistem koordinat Kartesian sebagai berikut. 


Inna) y (3.6) 

dengan mengingat bahwa x—rcos@ dan y-rsin0. Bentuk dari garis arus itu 
simetris terhadap sumbu y jadi kita hanya perlu mencari solusi untuk x 2 0 saja. Dari 
pertimbangan garis arus, cara terbaik untuk menghitung koordinat titik-titik garis 
arus adalah dengan menentukan nilai x, misalnya XS, dan menyelesaikan persamaan 
(3.6) untuk menghitung nilai y, untuk nilai x ditentukan, persamaan (3.6) dapat 


ditulis ulang menjadi bentuk f(y) — 0, dimana 


#(y) —y' —yy#(xs? —1)y 15lxs' ty)ml(xS? 4y?)- yxs? (3.7) 


Apabila nilai r diberikan, maka nilai fungsi fr) yang harus diberikan oleh persamaan 
(3.5) dapat dihitung. Demikian juga apabila nilai y diberikan maka nilai fty) dapat 
dihitung dengan rumus (3.7). Tetapi apabila nilai f yang diberikan, pertanyaannya 
adalah bagaimana cara menghitung r atau y ?. Untuk memudahkan pembahasan 
varuabel r dan y akan diganti oleh variabel bebas umum yaitu x. Masalah adalah 
apabila bentuk fungsi f(x) diketahui secara eksplisit, bagaimana cara menghitung nilai 
x seandainya f diberikan, sedangkan fungsi x(f) tidak diketahui bentuknya. Secara 
singkat kita akan membahas metoda numerik untuk menghitung solusi dari 
persamaan umum ftx)-0. Supaya dapat diselesaikan dengan baik, masalah harus 
dibatasi sebagai berikut. Bentuk f sebagai fungsi x diketahui secara eksplisit dan 
medan solusi dibatasi dari Xmin sampai Xmax yaitu Xmin Xx X Xmax. Ini juga 
berarti bahwa solusi XS pasti berada dijajaran harga 
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Xmin X XS X Xmax (3.8) 
Nilai Xmin dan Xmax ditentukan berdasarkan masalah yang dihadapi. Sebagai 
contoh batas medan solusi untuk persamaan (3.5) adalah Xmin-1 dan Xmax-3, 
walaupun kita dapat juga memilih Xmin-2.5 atau Xmax-4. Masalah dibatasi lebih 
lanjut dengan menuntut bahwa f(x) adalah suatu fungsi monoton, yaitu nilai f 
meningkat dengan meningkatnya nilai x (monoton naik) atau sebaliknya f selalu 
menurun dengan meningkatnya nilai x (monoton turun). Karena f(x) adalah fungsi 
monoton dan f(XS)-0, dimana solusi XS memenuhi persyaratan (3.8), maka 


f(Xmin).f(Xmax) «0 (3.9) 
Untuk fungsi monoton naik f(Xmin)» 0 dan f (Xmax)»0 sedangkan untuk 


fungsi monoton turun f (Xmin)» 0 dan f(X max) «0. Bagaimanapun juga persyaratan 
(3.9) selalu terpenuhi kecuali bila Xmin - XS, atau Xmax - XS, yang merupakan kasus 
remeh ("trivial") karena solusinya sudah diketahui. 

Pertama-tama ditentukan bahwa 


Xkiri — Xmin 
(3.10) 
Xkanan — Xmax 
dan ini berarti bahwa 
Xkiri « XS « Xkanan (3.11) 
f(Xkiri).f(Xkanan) « 0 (3.12) 


Selang (Xkiri, Xkanan) dimana solusi berada sekarang dibelah dua, yaitu kita menebak 
nilai solusi sebagai berikut 


Xtebak — Hxkiri 4#Xkanan) (3.13) 


Apabila f(Xtebak)-0 maka solusi sudah diperoleh yaitu XS-Xtebak,. Tetapi bila 
f(Xtebak) #0 maka kita harus melakukan pengujian untuk menentukan apakah XS 
berada diselang (Xkiri,Xtebak) atau (Xtebak,Xkanan). Apabila 
f (Xkiri).f (Xtebak) «0 ini berarti bahwa XS berada diselang (Xkiri,Xtebak) dan lebar 
selang solusi dimana solusi berada yaitu selang (Xkiri, Xkanan) dapat dipersempit dengan 
memperbaiki nilai Xkanan sebagai berikut 

Xkanan — Xtebak (3.14) 
Sebaliknya bila f(Xkiri).f(Xtebak) » 0, maka perbaikan nilai yang harus dilakukan adalah 
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Xkiri — Xtebak (3.15) 
Prosedur kemudian diulang lagi mulai dengan rumus (3.13) sampai akhirnya nilai solusi 
ditemukan secara eksak, yaitu f(XS)-0, atau bila lebar selang dimana solusi berada sudah 
bernilai cukup kecil sesuai persyaratan 

Xkanan- Xkiri « (3.16) 
dimana & adalah suatu bilangan positif kecil, yang kita tentukan sebagai ukuran ketepatan 
solusi, misalnya & —10“. Prosedur numerik di atas dikenal sebagai metoda pembelah 
selang atau metoda "bisection". Metoda ini selalu konvergen, tetapi membutuhkan jumlah 
iterasi yang cukup besar karena berderajat 1, yaitu galat hitungan dalam satu iterasi 
berbanding lurus dengan galat pada iterasi sebelumnya. Metoda ini sangat tidak efisien dan 
hanya dipakai oleh para pemula yang sedang belajar. Suatu metoda lain yang sedikit lebih 
efisien dan sangat mirip adalah metoda "regula falsi" atau metoda solusi palsu. Metoda ini 
pada dasarnya sama dengan metoda pembelah selang kecuali rumus (3.13) untuk 
menghitung solusi tebakan harus diganti oleh rumus berikut 


Kisbak 5 Kade 2 NAN ii (3.17) 
fkanan — fkiri 
dimana fkanan#f(Xkanan) dan fkiri-f(Xkiri). Kedua metoda di atas adalah contoh dari 
metoda dua titik atau metoda selang pengapit solusi. 
Suatu metoda lain yang jauh lebih efisien , yaitu metoda berderajat 2, adalah 


metoda Newton yang merupakan contoh dari metoda satu titik dan dapat dijelaskan 


sebagai berikut. 
Metoda Newton membutuhkan nilai tebakan awal yang disebut Xlama dan dapat dihitung 
sebagai berikut 

Xlama — 7 (Xmax- Xmin) (3.18) 


atau dengan rumus 


Sebebiaie Kiat 2 AKAN KAT iri 
fkanan — fkiri 








3 Xmax— Xmin 3 
Xlama — Xmin FKmax) Kemang 1 Ain) (3.19) 


Tebakan baru yang lebih baik kemudian dihitung dengan menggunakan rumus berikut 





f (Xlama) 
f'(Xlama) 





Xbaru — Xlama (3.20) 
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dimana f'(x) adalah turunan pertama dari fungsi f(x). Prosedur iterasi dengan metoda 


Newton dianggap telah konvergen apabila persyaratan berikut terpenuhi 


(Xbaru — Xlama| ce (3.21) 
apabila persyaratan (3.21) belum terpenuhi maka kita harus melakukan perbaikan nilai 
tebakan sebagai berikut 

Xlama — Xbaru (3.22) 


dan nilai Xbaru dihitung ulang lagi dengan rumus (3.20). Metoda Newton adalah metoda 
yang sangat efisien, tetapi tidak selalu berhasil atau konvergen, misalnya apabila nilai 
f(Xlama) dalam rumus (3.20) menjadi nol. Karena buku ini bukan buku tentang metoda 
numerik maka kita tidak akan membahas semua metoda yang ada ataupun membahas 
tentang masalah konvergensi dan lain sebagainya. Para pembaca yang ingin mempelajari 
masalah ini secara lebih mendalam dapat memperoleh informasi yang diinginkan dengan 
membaca buku teks misalnya yang ditulis oleh Hamming (1962), Carnahan Etal (1969), 
Gerald and Wheatly (1984), Chapra and Chanale (1985), Cheney and Kincaid (1985), 
Press Etal (1986), Nakamura (1991) dan lain sebagainya. 

Masalah numerik lain yang sering dihadapi adalah bagaimana cara menghitung 
integral, misalnya dalam menghitung nilai koefisien deret Laurent C(m). Penjabaran umum 
masalah ini dapat diberikan sebagai berikut. 

Data yang diberikan atau diketahui adalah himpunan titik atau pasangan nilai f dan x, 


misalnya (x(i),f(i),1 si IMAX). Titik-titik data (x,f) tersebut boleh jadi diperoleh 
sebagai hasil eksperimen atau dari suatu rumus analitis untuk f(x) yang tidak diketahui 
bentuknya. Dalam hitungan nilai koefisien deret Laurent, integral yang harus diketahui 
diberikan oleh rumus (2.18) yang pada intinya adalah 
Y(m) - 9ZA.exp(imO)d0 
Data yang diketahui adalah sejumlah koordinat titik ZA yaitu (XA(i),YA(i)) dan petanya 
di bidang-C yaitu @(i). Karena ZA adalah variabel kompleks, maka integral I(m) adalah 
satu bilangan kompleks yang dapat diuraikan menjadi komponen riil dan imajinernya, yaitu 
I(m) - ReI(m)--iIm(m) 


dimana 
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ReI(m) — $IxA- cos(m0) — YA.sin(m0) |.d0 

Im I(m) — $(xa.sin(mo) #-YA.cos(mO)|. d9 
Kedua rumus di atas dapat ditulis dalam bentuk umum sebagai berikut 

1- f(x).dx 
dimana x adalah 0, sedangkan f untuk komponen riil adalah 

f - XA.cos(m0)- YA.sin(m0) 
dan untuk komponen imajiner adalah 

f — XA.sin(m0) 4 YA.cos(m0) 
Karena jumlah koordinat aerofoil ZA yang diberikan itu terbatas, maka f yang diketahui 
juga terbatas. Sekarang marilah kita bahas kasus umum dimana kita ingin menghitung nilai 
integral dari fungsi f(x) dimana data yang diberikan adalah sejumlah terbatas pasangan 
nilai (x,f), misalnya (x(i),fG)) untuk i-1 sampai dengan i-IMAX. Untuk memudahkan 


pembahasan kita anggap bahwa selang Ax antara dua nilai x(i) yang berurutan itu bernilai 
tetap, yaitu untuk semua nilai i 


Axli) —xli#1)-xli)—Ax (konstanta). 
Misalnya kita ingin menghitung nilai integral dari f(x) antara x(i) dan x(it1) untuk 
ISIS IMAX -—1. Bila dalam selang nilai x tersebut kurva f(x) boleh dianggap dapat 
diwakili secara akurat oleh sebuah fungsi linier, yaitu 
Af @) 
Axti) 
maka integral dapat dihitung berikut 





f(x) fi) .x untuk xi) sx S x(i H1) 


il 


Uxi),x(i 41) - il 4 arti) .x “0)jax 


Axli) 





dimana Af (i) -f(i--1) -f(i) dan akhirnya 





ie -(O6 do) LAKI fa wr) 


atau 
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Ki,i--1) — Meta) #10). Ax6) 


Kurva f(x) diselang Axli) diberikan sebagai garis lurus, jadi rumus (3.23) sesungguhnya 
adalah luas trapesium dibawah garis f(x) tersebut dan oleh karena itu rumus (3.23) dikenal 


sebagai rumus trapesium. Rumus ini dapat diterapkan pada semua data yang diberikan dan 
dengan demikian integral untuk seluruh data adalah 


(li) --1)).ax(i) 


dan karena nilai Ax(i) tetap maka 


: 


I- 1, IMAX) -— 5 


U 


Axli) SU f(i-1) tri (@)) 


il 


1 


KAN (3.24) 


untuk kasus dimana f(1)-f(IMAX) rumus di atas dapat disederhanakan menjadi 


IMAX-—1 
I-Ax. Xf(i) (3.25) 
Rumus trapesium berbentuk sangat sederhana, tetapi mempunyai keakuratan hanya 
berderajat dua yaitu berbanding lurus dengan (Ax)”. Untuk mendapatkan rumus dengan 


derajat keakuratan yang lebih tinggi, yaitu derajat tiga atau sebanding (Ax), kita harus 


menggunakan polinom berderajat dua untuk mewakili bentuk fungsi f(x), yaitu 


fx) - A4B(x— xi )--celx— x)) 
Untuk menghitung koefisien A, B dan C kita butuh tiga titik data yaitu (x(i),f@)), 
(KG@#1),f(--1)) dan (x#2),f-2)). 

f(i)-A 

f(i-1) - AB.Ax tCAx? 

f(iw2)- A 42B.Ax -4CAx? 
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Dengan sedikit pengolahan aljabar akhirnya kita dapat memperoleh hasil berikut 
A-fli) 


—3f(i) -4f(i-1)-f(i--2) 
2Ax 








f(i)-2f(iw1)--f(i--2) 
2Ax 
Integral dari fungsi f(x) dari x(i) ke x(it2) adalah 


Cc 








it 


Kis) - Fasa x(i)) -clx— xD) Lax 


yang setelah diolah lebih lanjut akhirnya menjadi 


Wi,i--2) — Ani) war(i1)11612)) (3.26) 


Rumus di atas dikenal sebagai rumus Simpson 1/3. Perhatikan bahwa rumus tersebut 
berlaku untuk 2 selang x, dan apabila IMAX bernilai ganjil atau jumlah selang adalah 
genap maka 


Dai (2i—1) war(2i) -1(2i-4-1)) 


- V1, IMAX) — 


sal 


dimana IMAX2 — #(imax--1) 


Rumus integral di atas dapat ditulis ulang menjadi 


ia KG) sitaan) 15 (2112 Yio1-1) (3.27) 


i-l i-2 
Analisis galat yang lebih teliti dapat menunjukkan bahwa galat untuk rumus (3.26) atau 
(3.27) adalah berderajat 4 bukan derajat 3. Bila fungsi taksiran untuk f(x) adalah polinom 
derajat 3 yang melewati titik-titik i, it1, it2 dan it3 dapat ditunjukkan bahwa rumus 
untuk integral adalah 


KMAX2 - 5 (KMAX-1) (3.28) 


Rumus di atas dikenal sebagai rumus Simpson 3/8 dan mempunyai galat berderajat 4, 
seperti rumus Simpson 1/3. Tetapi perhatikan bahwa rumus ini berlaku untuk 3 selang jadi 
dapat digunakan untuk menghitung I(l,Imax) hanya bila (IMAX-1)/3 adalah bilangan 
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cacah. Seperti telah dinyatakan sebelumnya rumus Simpson 1/3 hanya dapat digunakan 
untuk menghitung I1,IMAX) apabila IMAX adalah bilangan ganjil. Seandainya IMAX 
adalah bilangan genap misalnya IMAX-2.JMAX, maka kita tidak dapat menggunakan 
rumus (3.27). Tetapi perhatikan bahwa bila IMAX adalah bilangan genap maka 
KMAX-IMAX-3 adalah bilangan ganjil dan nilai integral N(1,KMAX) dapat dihitung 
dengan Simpson 1/3 yaitu 


I1,KMAX) — 2.1) HF(KMAX)44 9" (21) 4 2 Oi v) 


Selanjutnya nilai integral IIKMAX,IMAX) dapat dihitung dengan menggunakan rumus 
Simpson 3/8 yaitu 


IIKMAX, IMAX) — SAK (KMAX) 4 3(KMAX 4 1)--3f(KMAX 4-2) «£ (IMAX)) 


Untuk IMAX bernilai genap, integral I1,IMAX) kemudian diberikan oleh rumus 
Y4, IMAX) — X1, KMAX) # II KMAX, IMAX) 


dimana KMAX — IMAX —3 dan KMAX2— T(KMAX-1) 


Rumus integral Simpson adalah rumus yang sering dipakai untuk menghitung suatu nilai 
integral secara numerik. Disamping rumus trapesium dan rumus Simpson sesungguhnya 
masih ada rumus-rumus numerik lainnya untuk menghitung integral, misalnya rumus- 
rumus Newton-Cotes, kwadratur Gauss, Integral Romberg dan lain sebagainya. Dalam 
buku ini kita tidak akan membahas rumus-rumus tersebut karena informasi mengenai hal 
itu dapat diperoleh dari buku-buku teks metoda numerik dan hal tersebut barada diluar 
batasan pembahasan atau isi buku ini. 

Masalah numerik terakhir dalam penerapan metoda deret Laurent yang masih 
perlu dibahas adalah interpolasi dan ini akan dibahas secara cukup rinci dalam bab berikut. 
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BAB IV 
INTERPOLASI UNTUK MEMPERBANYAK JUMLAH TITIK 
PADA PERMUKAAN AEROFOIL SEMBARANG 


Dalam masalah pemetaan konformal dengan rumus pemeta deret Laurent, data 
diberikan dalam bentuk sejumlah terbatas, JTA, koordinat titik-titik 
ZAAYAXAA),YAAG)) yang berada pada permukaan aerofoil ZA. Tetapi yang 
diinginkan adalah kemampuan untuk menghitung semua koordinat titik yang berada pada 
permukaan aerofoil ZA tersebut. Dengan perkataan lain kita ingin mempunyai suatu 
rumus analitis y(x) yang dapat diterapkan untuk menghitung ordinat titik, y, apabila nilai 
absisnya, x, diberikan. 

Pendekatan paling sederhana yang dapat dilakukan adalah dengan menganggap 
bahwa kurva y(x) disetiap selang x(j) antara XAA(j) dan XAA(j-H1) dapat ditaksir oleh 
fungsi linier 


Aylj) 
Axli) 


Misalkan absis dari suatu titik pada permukaan aerofoil ZA diberikan, yaitu XAI, 





y(x) — YAA(j) (x- XAA(1)) (4.1) 


dan nilai ordinatnya, YAI, harus dihitung. Pertama-tama kita harus mencari tahu diselang 
mana nilai XAI itu berada. Titik ZAA(1) dan ZAA(ITA) berhimpit dan berada di ekor 
aerofoil, sedangkan indeks j meningkat nilainya arah kebalikan putaran jarum jam. Ini 
berarti bahwa pada awalnya nilai XAA(j) mengecil dengan meningkatnya nilai indeks j 
sampai j-JHA dimana ZAA(JHA) adalah titik hidung (paling depan ) aerofoil ZA. 
Kemudian nilai XAA(j) meningkat dengan meningkatnya indeks j untuk j » JHA. Dari 
pengamatan di atas dapat diambil kesimpulan bahwa titik XAI itu berada dalam dua 
selang dengan nilai j yang berbeda. Selang yang pertama mewakili permukaan atas 
aerofoil sedangkan yang kedua mewakili permukaan bawahnya. Jadi untuk satu nilai XAI 
ada dua titik pada permukaan aerofoil yang memiliki XAI sebagai absisnya. Kenyataan ini 
menyulitkan penentuan diselang mana XAI berada. Untuk mengatasi masalah ini kita 
akan menggunakan pendekatan interpolasi parametrik sebagai rumus yang mewakili 
bentuk aerofoil ZA. Kita harus mencari parameter s yang berperilaku sedemikian rupa 
sehingga nilai s selalu meningkat dengan indeks j. Bentuk aerofoil kemudian diberikan 
oleh dua rumus interpolasi, yaitu x(s) untuk menghitung absis dan y(s) untuk menghitung 
ordinat titik yang diinginkan sebagai fungsi parameter s. Untuk satu nilai s kita menuntut 
agar nilai x(s) dan y(s) itu unik, yaitu hanya ada 1 nilai x dan 1 nilai y. Pertimbangan 


tersebut diatas akhirnya menentukan bahwa jarak atau panjang kurva dari titik ekor ke 
titik permukaan aerofoil yang lain sebaiknya dipilih sebagai parameter s. 

Apabila bentuk aerofoil diwakili oleh polinom linier sepotong-sepotong (untuk 
setiap selang) maka panjang kurva SAAG) dari titik ZAA(1) ke ZAA(j) dapat dihitung 
sebagai berikut : 





SAA() —0 
SAA(j) 5 SAA(j -—1) #AS(j-1) untuk 2S js JTA 
AS(j-1) - /Ax(j -1) 4 Ayli - 1) (4.2) 


Axlj —1) — XAA(j) — XAA(j—1) 
Ay(j-1)- YAA(j)- YAA(j 1) 





Rumus diatas dapat diterapkan untuk menghitung semua nilai SAA(j) sehingga 
sekarang kita memiliki himpunan harga (XAAG),SAAG)) dan (YAAG),SAAG)) untuk 
1S js JTA dimana nilai XAA(j) meningkat dengan indeks j. 

Letak titik yang ingin dicari koordinatnya, yaitu ZAI, ditentukan oleh panjang 
kurva dari ekor ZAA(1) ke ZAI sepanjang permukaan aerofoil, yaitu SAI. Nilai SAI 
diberikan sebagai data. Untuk mencari koordinat ZAI pertama-tama kita harus 
menetukan diselang s(j) mana nilai SAI berada. Karena s adalah fungsi monoton naik dari 
indeks j, maka prosedur yang diperlukan menjadi sangat sederhana. Dimulai dari j-1 kita 
uji apakah SAI memenuhi persyaratan berikut : 

SAA(j) S SAI S SAA(j —1) (43) 

Apabila SAI tidak memenuhi persyaratan di atas maka indeks j dinaikkan nilainya 
menjadi j-2, j-3 dan seterusnya sampai akhirnya persyaratan (4.3) terpenuhi. Misalkan 
nilai j pada saat itu adalah JI. Absis dan ordinat ZAI kemudian dapat dihitung sebagai 
berikut : 
Ax( JI) 
XAI — XAA (JI) « —“ (SAI—SAA/JI)) (44) 


ASI JI) 


Ay(JI) 
AS( JI) 


Prosedur interpolasi di atas memang dapat dilakukan dengan mudah, tetapi hasilnya tidak 


YAI — YAA (JI) 4 (SAI —SAAJI)) (4.5) 
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terlalu memuaskan karena kurang akurat. Permukaan aerofoil sesungguhnya berbentuk 
suatu kurva lengkung yang mulus disemua selang, tetapi ditaksir oleh kurva linier 
sepotong-sepotong atau oleh sepotong garis lurus di masing-masing selang. Hasil yang 
lebih akurat dapat diperoleh dengan menggunakan kurva kwadrat atau polinom 
berderajat lebih tinggi sepotong-sepotong. 

Misalkan diselang j di permukaan aerofoil ditaksir oleh sebuah polinom kwadrat. 
Polinom tersebut harus melewati titik-titik ZAAG) dan ZAAG-#). Tetapi polinom 
kwadrat melibatkan 3 koefisien yang harus dihitung nilainya dan ini membutuhkan satu 
data tambahan, misalnya titik data ZAA-1) atau titik ZAAGt2). Bentuk umum rumus 
polinom kwadrat adalah : 


y-A-#Blx— XAA(j)) clx — xAA ()) (4.6) 


Karena polinom ini melewati titik-titik j, j-1, dan jt2 maka : 


YAA(j) -A 
YAA(j--1) -A-4B.Ax(j) # C.Ax( 7) 
YAA(j 42) - A4B.(Ax(j) # Ax(j1)) 4 C.( Ax(j) 4 Ax(j 1) 


Pengolahan aljabar selanjutnya pada akhirnya dapat memberikan nilai koefisien-koefisien 
A, B, dan C sebagai berikut : 





A- YAA(i) 
. Ayli)tAylit1) | Ayljt1) Ayi) 
BA) Anji) (aa 3) id 





Ay(jt1) Ayli) 
Axti#1) Axti) 
Axlj)-Axlj-t1) 

Panjang kwadrat diselang j dapat dihitung sebagai berikut : 








jd 


jl 
AS(j) — Jas — Nax 4dy 
j i 


atau 
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AS(j) — | (aja 


x 
Turunan dari persamaan (4.6) adalah : 


SY 2 420x—xaA(i) 


jadi 
PA 
AS(j) - ik #(B420(x—xAA))))) | dx 
i 
Apabila ditentukan bahwa : 


f(x) (into x20ta- xaA GT 


maka integral diatas dapat dihitung dengan rumus Simpson 1/3 sebagai berikut : 
FN! 5 , ai : 
As) pax (G) sar is) on) (4.8) 


(i) —£(xaafi)) - (wap? 


14 
li "1 - 1(xaa() 0) - (uts #c.Ax(i)Y | . 
, : na 

f(i1) —1(xAA(j -1)) —fir(ms20.Axdi)) | 
Setelah AS) dihitung, kemudian SAAG) dihitung dengan rumus : 

SAA(j--1) -SAA(j)-AS(j) untuk 1 £ js JTA -1 (4.9) 

SAA(1)-0 
Perhatikan bahwa rumus (4.7) itu berlaku untuk 1« js JTA-2 dan untuk 
JHA « j X JTA -2 dimana JHA adalah nilai indeks j untuk titik hidung aerofoil. Untuk 
JFJHA-1 dan j-JTA-1 kita harus menggunakan polinom yang melewati titik-titik j-1, j 
dan j1, sehingga rumus untuk koefisien-koefisien A, B dan C harus diubah menjadi : 


A- YAA()) 
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Axlj 1) .Ay(i) « Axlj) .Ay(j—1) 
Axl). Ax(j—1)-(Ax(j —1) 4 Ax(j)) 


| Ayli) 23) 





B- (4.10) 





Axlj) Axlj-1) 

Axli—1)-t Ax(i) 
Himpunan data (XAAG),SAA()) dan (YAA()),SAA()) kemudian digunakan untuk 
menentukan rumus-rumus interpolasi x(s) dan y(s) yang mewakili bentuk aerofoil ZA. 





Tentu saja untuk kasus ini kita harus menggunakan rumus interpolasi polinom kwadrat 
sepotong-sepotong, yaitu untuk s diselang j, 


x(s) — XAA(j) #D(s—SAA(j)) -E(s— SAA(j)) (4.11) 


y(s) — YAA) #F(s - SAA(j)) - G(s- SAA(|)) (4.12) 
dimana D, E, F dan G adalah : 


D-— 





(4.13) 


Axli) # Ax(j 1) | Axlj #1) 23) 
AS(j)--AS(j-1) (AS(j-1) AS(j) 








aa H1) 29) 
£ AS(j--1) AS(j) 
AS(j) - AS(j--1) 





(4.14) 








)--ay( 1) ea “| (4.15) 





AS(j--1) As(j) 

AS(j)--AS(j-1) 
Rumus diatas berlaku untuk 1X j X JTA -2, sedangkan khusus untuk selang j-JTA-1 
nilai D, E, F dan G adalah : 


aa 3) 





(4.16) 





AS(j —1)'.Ax(j) - AS(j)”.Ax(j —1) 
13 ASASI “Mas -1 rAS()) an 
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AS(j —1)--AS(j) (4.18) 

AS(j —1)'.Ay(j) - AS(j) -Ay(j-—1) 
F-A8(ihas(i AS —1) #ASG)) (4.19) 
.Iasti) Asti-1) ta 





AS(j —1)  AS(j) 
Rumus interpolasi derajat 3 atau lebih tinggi dapat diturunkan dengan pendekatan yang 
sama walaupun rumus yang terlibat akan berbentuk jauh lebih rumit. 

Rumus interpolasi yang telah dibahas adalah beberapa contoh sederhana dari 
prosedur interpolasi pada umumnya. Menurut Carnahan etal (1969) fungsi interpolasi 
yang memenuhi persyaratan sebagai fungsi penaksir yang baik adalah fungsi berbentuk 
deret dimana suku-sukunya adalah anggota salah satu dari tiga himpunan fungsi 
sederhana yaitu : 

1. Himpunan fungsi monomial £(x9), n - 1,N) 

2. Himpunan fungsi tridiagonal ((cos nx, sin nx), n—1,N) 

3. Himpunan fungsi eksponensial £(e1X,nX), n—1,N) 

Dibawah ini kita akan membahas interpolasi berbentuk polinom, yaitu deret dengan 
monomial sebagai sukunya, secara agak lebih terperinci. 

Fungsi interpolasi yang baik harus bersifat mulus dan berbentuk analitis sederhana 
sehingga rumus untuk turunan dan integralnya dapat diturunkan dengan mudah. Fungsi 
interpolasi juga harus mampu menaksir harga fungsi yang sebenarnya dengan akurat dan 
harga tersebut harus dapat dihitung dengan mudah dan efisien, Persyaratan keakuratan 
dapat dipenuhi bila fungsi interpolasi itu paling tidaknya dipaksa harus melewati semua 
titik data yang diketahui, misalnya ((X55Y1), ILN). Bila fungsi f(x) adalah fungsi 
interpolasi yang dipilih, maka persyaratan tersebut diatas dapat ditulis sebagai berikut : 


f(x) -y(x) | 


(4.21) 
f(x), untuk i-1,2,...,N 


Polinom berderajat m dalam variabel x dapat ditulis sebagai : 
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B,(x) — Ya. (4.22) 


atau 


P (x)—a, ta, .xta,.x4..ba ok" 4.23 
m 1 2 3 m-1 


Polinom (4.23) bersifat mulus dan mempunyai turunan ke satu sampai ke-m, yang 
semuanya bersifat mulus, dengan rumus sebagai berikut : 





d" ( ) & stel (i sz 1)! xi 


—P,(x : ... (4.24) 
dx na (j-1-n)! : 
Rumus untuk integralnya adalah : 
mi a 
J P,(x)dx — konstanta 4)” (2) x (4.25) 
AaANJ 
Harga P,n(x) untuk x—X dapat dihitung dengan rumus (4.23) yaitu : 
P,(X)—a, ta,X ta,X4-...ta,,,X” (4.26) 


, 3 Ka aan : 1 : 
Cara hitungan diatas itu tidak efisien karena membutuhkan yet FI) perkalian dan 


dapat diperbaiki dengan menggunakan cara pengurungan berlipat atau cara Horner yang 
hanya membutuhkan m perkalian saja, yaitu : 


m 


B,(x) — (ay #X.(a, #X. (ay. 4X. (any #Xua, #X.a..,1)).-))) 427 
atau 
Pa (1) 5 (((-e (aya Xta,,) X-ray, 1). Xa-.ta,).Xa, X2) 420 


Polinom yang memenuhi syarat (4.21) adalah polinom berderajat (N-1) dan koefisien- 
koefisiennya memenuhi persyaratan berikut : 


N 
Hak, untuk i - 1,2, 3,...N (4.29) 
ja 
atau ditulis dalam bentuk perkalian matriks : 
(cu klas -Pxl untuk i,j 1, 2,3, ..N (4.30) 
dimana 
Cc, - xXx 
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Persamaan matriks (4.30) dapat diselesaikan untuk menghitung nilai aj, tetapi proses ini 
tidak efisien dan cenderung melibatkan galat atau kesalahan numerik dalam nilai koefisien 
aj. Suatu cara yang lebih baik untuk menghitung nilai a, dapat dijabarkan sebagai berikut. 

Bila polinom P,(x) dibagi oleh (x-X:) maka hasilnya adalah sebuah polinom 
berderajat 1 lebih rendah daripada m, yang dikenal sebagai polinom susut, dan suatu 
bilangan Sj yang merupakan sisa dari pembagian yaitu : 

P.(x)-P, , (x).(x— X,) 4S, (431) 
Persamaan (4.31) menunjukkan bahwa Sj sesungguhnya adalah nilai polinom P, (x) 
untuk x—Xj. Kenyataan ini dapat dimanfaatkan sebagai berikut : 

Pertama-tama polinom P, (x) dibagi Xj dan hasilnya adalah : 


PB, (x) —P, ,(x).(x— X,) 4S, (4.32) 
Karena P, (X,) - Y,, maka : 

SY (4.33) 
Selanjutnya polinom P,, ,(x,) dibagi oleh (x-xX,) dengan hasil 

Pr.i(x) “Pr. (x)-(x—X,) 5, (4.34) 
Pemasukan (4.34) ke (4.32) memberikan hasil berikut : 

B, (x).(P, »(x)-(x— x,)4-S,).(x—x,) 5, (4.35) 


Mengingat bahwa P, (X,) - Y,, maka persamaan (4.35) dapat disederhanakan menjadi : 
X, 58,.(X, — X,)45, 


atau 
Y,—$ 
S5, -—-——— 4.36 
Aa REA (4.39 

Pembagian P, ,(x,) oleh (x- X,) memberikan hasil berikut 

P,.2(x) - Br 3(x).(x— X3) 45, (4.37) 
yang bila dimasukkan ke (4.34) memberikan hasil 

P,..(x)-(P,..(x)-(x si x,) tS,).(x - x,) 4S, (4.38) 
Pemasukan (4.38) ke (4.32) akhirnya memberikan 

B, (3) -((P, (1). (x— X,) 55). (x— X,),).(x—X,)-t5, (439) 


Karena P, (X4) — Y,, maka 
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(S3: (X3 —X,) #8, ) (X5 — X,) 4S, 5 Y, 


atau 


3) 
IM 
S5 — UU (4.40) 
(Xx —X,) 
Prosedur diatas dapat dilakukan lebih lanjut dengan pembagian oleh 
(x- X4) (re) aa, (Xu) dengan hasil berikut : 





(23 -s, 
XX, 
2) 2-8, | (rx 4.41 
S, aa 3 (Ke -X,) (441) 
(aa -s, 
XX 
S. — OK) S, | /xe-x,) I-S, | /(xe-x.) | (4.42) 


dan pada umumnya untuk m 22 


Ora) 
| m1 1 -s,| 
Xi Xi Sa (re 


mat “|| se (Xx 





Kn X) (4.43) 


mil” Am 
ml 


Polinom Pyy(x) dapat ditulis dalam bentuk menggunakan koefisien-koefisien Sj untuk 
1S ix m-#1 sebagai berikut : 


P, (x) -(-x((Snsae(x— Xp), )-(x— X4) #8, ):-)lx—x,) 5, (444) 


Seringkali disamping menginginkan harga f(x) yang ditaksir oleh Pyn(x) kita juga ingin 
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menghitung nilai f'(x) yang ditaksir oleh turunan P,a(x) yaitu Te). Nilai turunan ini 


juga dapat dihitung dengan mudah sebagai berikut : 
Pertama-tama kita rumuskan fungsi Tyy(x) 


Tx) Si: (x— X,,) 5, (4.45) 
Tx) Tx): (x— X1) ta. (4.46) 
Tx) To.ix).(x— X,) 4S, (4.47) 
Perumusan diatas menunjukan bahwa Tym(x)“Pm(x) dan turunan Ty,y(x) itu sama dengan 
Pm): 
Sekarang perhatikan bahwa : 
d d 
—T) - Tenar) #(x— X,)-—T,, ,(x) (4.48) 
dx dx 
d d 
Tor) Tax) #(x—X,)- Tax) (4.9) 
dx dx 
d d 
— Tx) - Tx) #(x— Xu, 1) -— Tx) (4.50) 
dx dx 
TR) Sa. (4.51) 
dx 
Bila ditentukan bahwa untuk 1X js m 
V3) - LM) (4.52) 
j dx 1 : 


maka turunan-turunan Tx) dapat dihitung secara beruntun sebagai berikut : 


Ur (Xx) —Senra (4.53) 
U,(x) — Ti(x) #(x—X., 4)-U, (x) (4.54) 
Uzlx) — T,(x) #(x—X,, ,)-U,(x) (4.55) 
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U,ra(x) 5 Tr.2(x)t(x— X,).U,.o(x) (4.56) 
U,, (x) 5 T.y(x)-#(x— X,)-U,,.i(x) (4.57) 


Karena Tyn(x) - Pm(x) maka 

d 

—P,(x)-— U,,(x) (4.58) 

dx 
dan seperti dapat dilihat nilai U,,(x) dapat dihitung dengan rumus rekursif atau rumus 
beruntun (4.53) sampai dengan (4.57). Perhatikan bahwa nilai Ty5(x) dapat dihitung 
secara beruntun dengan rumus (4.45) sampai dengan rumus (4.47). Untuk kasus selang 
seragam dimana untuk semua j 

Ax 5 Xi, —X, 5Ax (konstanta) 


rumus-rumus yang dijelaskan diatas dapat disederhanakan sebagai berikut : 


SY, (4.59) 
s,- An (4.60) 
s, Baen (4.61) 
8, Jne Ata Ay) (4.62) 


s3 1 (Ay, —3Ay, 43Ay, —Ay,) 
24 Ax' 
dan seterusnya. 





S, (4.63) 


Dalam rumus (4.60) sampai dengan (4.63) dirumuskan bahwa : 

Ay, Sia Ih (4.64) 
Sekarang marilah kita rumuskan beda hingga Newton maju dan mundur sebagai berikut : 
Beda maju: 

Ny, 5y, untuk j-1, 2, ..., m, mt1 


(Aya -Ny,) 


Ny, - untuk j-1, 2, ..., m 
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(Aya -A'y,) 


Ny, - 5 


untuk j-—1, 2, ..., m-1 


j2 Xp 


(Aya SA) 


Ay, - untuk j-1, 2, ..., m4t1-k (4.65) 
(rom ag) 


KA ng aan 

Ay, — | Ysa y)) untuk j-1 
Xiam “Xj 

Beda mundur : 


V5, untuk j-—1l, 2, ..., m, mt1 


(vsy, 2 V'y) 


Vy,- untuk j- 2, ..., m#1 
Xj Xp 
Viy. vi 
V'y, Ah y-Vya) untuk j-3, ..., mt1 
Xj Xin 
Vy yi, 
Vy, | Yi Ya) untuk j-k41, k42, ..., mt1 (4.66) 
Xx X4 
v4 23 yud Y 
v"y, ( Yi Ys ) untuk j-m-1 
Xp Xom 


Rumus (4.23) dan (4.36) sekarang dapat ditulis dalam notasi beda maju sebagai berikut : 
S5 Ay, (4.67) 
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WN NY, —AY 
ME KK 
Rumus (4.40) dapat diolah lebih lanjut menjadi : 





S, 








AY, (4.68) 























S, XX Ay, 
(X, Xx, CX, x,) X3 Xx, 
s5) N) Ay) 
3 
(X, ai XX, rw x,) 
gi 1 (as winda KI Ad | 
: (XX) i X3 —X, 
1 
Ka ia 
3 ag Y? y) 
S, Ay, (4.69) 
Pada umumnya dapat ditunjukkan bahwa 
Sari 5 A Yura (4.70) 
Rumus untuk S dapat juga ditulis dalam notasi beda mundur, yaitu 
Sei SV Yura (4.71) 


Fungsi polinom (4.44) dapat ditulis dalam notasi beda maju dan beda mundur sebagai 
berikut : 
P, (x) — (-.Laryulx —X,) HAM, Jk X4) 4 Ay, hAty ix -X,) Ay, 

(4.72) 


P, (x) — - : (rasa x — Xi  Yy Jk X1) Las 22 £ #Vyh (x- XV, 
(4-73) 

Polinom dalam bentuk (4.72) dikenal sebagai polinom beda maju Newton.sedangkan 

(4.73) dikenal sebagai polinom beda mundur Newton. 

Fungsi interpolasi polinom dapat juga diberikan dalam bentuk polinom Lagrange, yaitu : 

Polinom Lagrange derajat 2 
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(x- Xxx.) (x-X)(x—x,) 


NY, (as Tila xXx) Ng 
(4.74) 


ae An 


2 " (Xp — X,) 








yang dapat ditulis ringkas seperti berikut : 


P,(x) — Ia P3 x) NA use x,) Ny pe) NX, (415) 


ui ta “2 1 (aa— “3 - —X 
atau 
Nb Ja da x) Xx (4.76) 


i- Ta (&-x,) 


Polinom Lagrange derajat m dapat ditulis sebagai berikut : 


wal aa (x—x)) 

P W- Tx) x (4.77) 
Polinom Newton dan polinom Lagrange sesungguhnya adalah fungsi yang sama, hanya 
ditulis dalam bentuk yang berbeda. Bentuk polinom Lagrange itu cukup rumit dan nilai 
Pm(x) tidak dapat dihitung secara efisien. Nilai turunan Py,(x) juga sangat sulit untuk 
dapat dihitung secara efisien. 

Sekarang marilah kita terapkan hasil-hasil di atas untuk memperbanyak jumlah 
titik pada permukaan aerofoil ZA. Pertama-tama kita harus menghitung distribusi 
panjang kurva SAA(j) seperti telah dilakukan sebelumnya dengan rumus (4.2) untuk 
taksiran linier atau rumus (4.6) untuk taksiran kwadrat. Di sini kita akan menggunakan 
polinom Newton derajat 3 yang harus melewati 4 titik data. Untuk kurva diselang j 
dimana 2 £ j X JHA-2 dan JHAt1 SXj X JTA-2, polinom penaksir harus melewati 
titik-titik j-1, j, j#1, Jt2. Rumus interpolasi untuk kasus ini adalah : 


P, (x (x ja IK Yi- lx Ka) t Aya lx 3 x) t Aya (x sei x,) t Ayi (4.78) 
Seperti telah dijelaskan sebelumnya (lihat rumus (4.8)) yang ingin dihitung nilai 
sebenarnya adalah turunan polinom dititik j, j#-1/2 dan j--1. 
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Untuk menghitung nilai turunan-turunan tersebut kita perlu merumuskan 
beberapa fungsi pembantu sebagai berikut (lihat persamaan (4.45) sampai dengan (4.57)). 


Tx) 5 Ayla 3 5) A Ryan 
Tx) - Tx)x 3g! x) 23 Aa 


T,(x) — T,(x)(x - Te) # 3 


U,(x) — an) Ayi 
U,(x) (3) — Ula x,) 4 (3) 





U, (x) - F5) — U,(x)(x Ke) #T(x) 


Karena T,(x) — P,(x) maka rumus untuk turunan yang diinginkan diberikan oleh U,(x), 


yang dapat dihitung secara beruntun dengan bantuan rumus-rumus di atas. Untuk tujuan 
pembuatan program komputer rumus-rumus di atas dapat diringkas sedikit menjadi : 


aloe ae X) 4 Aya 
x)(x - x) #AYi, 
U, (x) — Aya lx — x) 4T,(x) 
U,(x) — U, (x)(x — Ss) 4 T,(x) 


— 
Ex 

wd 

11) 
b, 

LN 


— T, “ 
(4.79) 


Untuk interpolasi diselang j-1 dan j-JHA, polinom penaksir harus melewati titik-titik j, 
JL, jt2, dan jt3 sehingga rumus untuk polinom adalah 


P, (x) — IKA — Kesal 4 #y, |(x - Xu) 4 Ay, (x — x) 4, (4.80) 


Untuk kasus ini nilai turunan dapat dihitung dengan rumus-rumus rekursif berikut : 


Tx) - Ayy(x EF Xysa) tAy, 

Tp3) 5 (xxx, ) HA, 
U,(x) Ay, lx X1) # Tx) 
Us(x) —U, (x)(x— x,) 41, (x) 

Untuk selang j-JHA-1 dan j-JTA-1 rumus untuk polinom penaksir adalah : 
P,(x) (Ket —X) Sy, ae Xt Ay, x 2x aja 682) 


karena kurva harus melewati titik-titik j-2, j-1, j dan j#1. 


(4.81) 


Rumus-rumus rekursif untuk menghitung turunan dalam kasus ini adalah : 


- Ayat X) Ay, 
- a)la-X,,) Ay,» 
U,(x) Ay, (x—X,,) To (x) 
Uglx) — U, (x)(x— x,.,) 1, (x) 


(4.83) 


Setelah nilai-nilai turunan dihitung, kemudian distribusi panjang kurva SAA(j) dihitung 
dengan menggunakan rumus (4.8) dan rumus (4.9) dengan catatan bahwa nilai f(x) 
diberikan rumus 
1 
19 -it us) (4.84) 

Selanjutnya rumus-rumus interpolasi bentuk aerofoil, yaitu x(s) dan y(s) dapat dihitung 
dengan menggunakan himpunan titik data (XAA(j), SAA(j)) dan (YAA)),SAA)). 
Bila rumus interpolasi yang digunakan adalah polinom Newton maka rumus-rumus yang 
dibutuhkan adalah : 
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Untuk j-1 


Ts) — A'x.(s -$,,,) 4x, 
Ts)  T,(5)(s—S,,,)-“A'x, (4.85) 
x(s)— P,(s) — T,(5)(s -s,) 4X 


T,(s) - A'y,(s—5,..) 4Ay, 
5) T(5(5—S,,,) Ay, (4.86) 


y(5) -P(5)  T5(5)(s—S,,,) 4, 
Untuk 2 jx JTA-2 


T,(s) - Axya(s—S,..) Au, 
T,(s) - T(5)(s—S,) HA'x, , (4.87) 
x(s) — P,(s) 5 T,(5)(s ol at Xia 


DS) Ayyals Sy.) Ay, 
Tas) - Ti (5)(s—S,)--A'y,., (4.88) 
Y(5) Pas) - T5)(5—S,.,) 4, 

Untuk j — JTA-1 


Ts) -— Nx,a(s - s,) tAx, 


Bs) -T5)(S—S,.,) 4A'x,., (4.89) 
x(s) — P, (s) 2 15)( “aj tX, 
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Ts) - Ny,as 3 s) 21 Ry,» 
T,(s) - 5)(s-S,,) HA, (4.90) 


y(s) — P,(s) — T(5)(s s0 tY 

Suatu kelemahan dalam interpolasi dengan fungsi polinom seperti telah dijelaskan di atas 
adalah kenyataan bahwa nilai turunan di setiap titik data j bila didekati dari kiri itu 
berbeda dengan nilainya bila didekati dari kanan. Dengan perkataan lain distribusi 
turunan tidak mulus dan patah-patah di setiap selang. Kelemahan ini dapat diatasi dengan 
menggunakan polinom lentur ("spline") dimana kita menuntut bahwa turunan fungsi 
harus bersifat mulus. Polinom yang akan dibahas dipilih berderajat 3 dan harus memenuhi 
persyaratan berikut untuk setiap selang j. 

1. Polinom harus melewati titik-titik ujung selang j dan j-1. 

2, Disetiap ujung selang turunan pertama dan kedua dari polinom penaksir 

harus bersifat mulus. 

Karena turunan kedua polinom lentur harus bersifat mulus maka turunan kedua ini paling 
tidaknya harus berbentuk fungsi linier sepotong-sepotong. Untuk selang j rumus untuk 
turunan kedua adalah : 


" hai " 
f 0 -—-e-—x)at (4.91) 


jl j 


Persamaan di atas dapat diintegralkan 2 kali berturut-turut dengan hasil berikut : 


1 f3 f3 





f(x) Sa kl tax) ex, (4.22) 
f(x) - : La (ex) 4 tita xXx) Kila—x)aK, (4.93) 


Untuk x“Xj nilai turunan dan fungsi adalah fj' dan fj, sedangkan menurut persamaan 
(4.92) dan (4.93) nilai turunan dan fungsi tersebut adalah K, dan K,. Ini berarti bahwa 


Kar 
K,-f, 
dan rumus (4.92) dan (4.93) dapat ditulis menjadi 
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| f3 —ih 2 
Po Sa) Ala) 4.94 
Ox ara) (434 
Meh f3 Sea, 2 
f ja 3 Xx f3 Jaa), 4.95 
(x) 23 (x 1 1G Fa x) t 3(x J4 i (4.95) 


Untuk ujung kanan selang j, yaitu titik j-1, persamaan (4.95) memberikan hasil berikut : 


1 1 
Nai (ita (ada Jang tan, 








dimana x, “X1 -X, 


Rumus di atas dapat disederhanakan menjadi 


fe. -f, 


i - Hg 


Ax Ga 21 )Ax, Hj (4.96) 
j 


untuk selang j-1, rumus interpolasi f(x) berbentuk mirip persamaan (4.95) kecuali bahwa 
indeks j dalam (4.95) perlu diganti oleh indeks j-1. Ujung kanan selang j-1 adalah titik j 
dan rumus (4.96) dengan mengubah j menjadi j-1, yaitu : 


ff, 1 
F- 1 — la 4211) Ang H3, (4.97) 


Pengurangan persamaan (4.96) dan (4.97) dapat memberikan hasil berikut : 





fi —f 


AI 


-f, 
AD 493) 
Ary, 


Sekarang perhatikan bahwa rumus untuk turunan (4.94) bila diterapkan pada selang j-1 


(5 ta Jan, (sar) 


menjadi : 


Fa) Maia 


2 eh An laen najar, (4.99) 
j- 


Khusus untuk ujung kanan selang j-1, yaitu titik j, rumus (4.99) memberikan hasil berikut 
1 
f4 2 al (Jan. HAN, 
yang dapat disederhanakan menjadi : 
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ff Haa H3) An, (4.100) 


Penggabungan rumus (4.98) dan (4.100) akhirnya dapat memberikan persamaan 
dibawah ini : 

ff, £-f, 

Ap HJ An 4 Hi H2 )Ax, 56) Sa 

(rama, sala sam)an 2g an Ab 


yang dapat ditulis ulang menjadi : 


ff ff 
Ax, f1 42 Ax,, An Ji" 4 Ant, —6 | Ld 
Ax Ax 


jatja | (4.101) 
i j4 

Persamaan kemulusan turunan kedua (4.91) dapat juga diberikan oleh persamaan berikut: 

f3 f" 

f(x) -— Saya —x)ari, (4.102) 


yang dapat diintegralkan dua kali dengan hasil 


j 





f(x) — Lai (x x) fu (x x) Ai (4.103) 
“jua ja ja ja . 
2 Ax 
Lt fi 1 5 
f(x) ——— 3 dapa x) 4 mal, —3) Sa Saat, (4.104) 





6 Ay 
Untuk x — X, persamaan (4.104) dapat ditulis menjadi 


LE -- Har HA! Ax, (4.105) 


ja “—N— jua 


sedangkan untuk selang (j-1) dititik X, , diperoleh hasil berikut 

U (op tt ) Mx, 41. Ax 4.106 
Sia 5 H3 Ang, H.A, (4.106) 
selisih (4.105) dikurangi (4.106) dapat ditulis sebagai berikut 


Ii Iis (4.107) 


fl —f — H2 Pe) An, (or ag) Axy Jt AT 
j j- 


sedangkan persamaan (4.103) untuk x — X, dapat ditulis menjadi bentuk berikut 
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1 
Fira — j 


f'- ata H3). Ax, (4.108) 


jl 


selisih dari persamaan (4.107) dan (4.108) dapat ditulis menjadi bentuk berikut untuk 
25X ITA-—1 


ha K-f. 


Ang yet 4 AA, , # An)" 4 An, — ga Maa 


Sal (4.109) 


Perhatikan bahwa persamaan (4.109) itu tepat sama dengan persamaan (4.101) 
karena keduanya memiliki persyaratan kesinambungan turunan kedua dari fungsi 
interpolasi, walaupun diturunkan dengan pendekatan yang sedikit berbeda. Sistem 
persamaan (4.109) melibatkan JTA variabel anu yang harus dihitung nilainya, yaitu f" 
untuk IS jXJTA tetapi hanya terdiri dari JTA-2 persamaan. Oleh karena itu 
dibutuhkan dua persamaan tambahan untuk menyelesaikan persamaan tersebut secara 
simultan. Kedua persamaan tambahan yang dibutuhkan pada dasarnya adalah informasi 
tentang nilai f," dan ft, yang tentu saja tergantung pada masalah yang dihadapi. Secara 
umum dapat dikatakan bahwa ada empat jenis masalah kondisi batas yang berbeda, yaitu 
Kasus I : 

Untuk kasus fungsi kubik lentur alamiah turunan kedua fungsi interpolasi di kedua 
ujungnya bernilai nol, yaitu : | 

fa, 50 (4.110) 
Kasus II : 
Kasus yang bersifat cukup sederhana adalah dimana nilai turunan pertama di kedua ujung 
bernilai nol, yaitu : 

fi 50 (4.111) 

Penerapan persamaan (4.96) di titik j-1 memberikan hasil berikut 


ct - - aa 
Ta Ax, 


untuk f," — 0 persamaan di atas dapat diolah menjadi 


Da P5 


Kh f 
f — 1 pu tah 
2 Ax? 


Penerapan persamaan (4.105) di titik j-JTA-1 memberikan hasil berikut 


(4.112) 
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u 1" frra 25 frra | LN 
Ara Haa 5 6, Hb 
JTA-—Y X yrA 





Untuk fl, — 0 persamaan di atas dapat diolah menjadi 


1 ni 
ht -2 " 2 
fira pora 3 sa 1 (4.113) 





Kasus III : 
Untuk kasus umum dimana fungsi yang ditaksir tidak bersifat periodik dan nilai f" serta 
Fira tidak diketahui, kita harus membuat taksiran sebagai berikut. 
Ujung kurva dianggap dapat ditaksir cukup baik oleh suatu polinom derajat 3 
yang melewati 4 titik data. Polinom Newton beda maju derajat 3 yang melewati titik j-1, 
2, 3 dan 4 adalah 
p(x) (A5, (x—x,) 5). (x—x,) A15 Ya), 
dan turunan pertamanya di titik j-1 adalah 
5 (Af (Ax, & Ax,) — A4, ).Ax, HA, — K, (4.114) 
Polinom Newton beda mundur derajat 3 yang melewati titik-titik j-JTA-3, JTA-2, JTA-1 
dan JTA adalah 
P(x) “(Va (x— Xgp-2) # Von): (x— Xora-1) Vira N x— Kura) #-Hara 
dan turunan pertamanya di titik j-JTA adalah 
Pora Vga (Vxyra 1 # Vxyra ) # Vforr ) Vxyr 1 #Vifry 5K, (4.115) 
Arti dari simbol beda hingga adalah sebagai berikut 
An, - 2 2 
At, (to,u —$,)/Ax, 
Af, 5 (Af, — At, )/(Ax, #Axy,,) 


tAx,:) 


jl 


(a 

3 5 — N1,)/(Ax, 4 Ax 
pr -t)/Ax., 

r (Vir — Vir, ,)/(ax,., # Ax,.)7 

(ver, — var,)/(Ax,.. # Ary, Ary.) 


.5 


51 


setelah f dan fi, dihitung nilainya, yaitu K1 dan K2, kemudian persamaan (4.96) dan 
(4.105) diterapkan sebagai berikut 











pi Lgu ah tas (4.116) 
2 Ax, Ax, 
1 fyra —£, K 
fa fu 3-3 “MA 3 2 4117 
JTA 2 JTA-1 Ra Xena ( ) 


Kasus IV : 
Dalam kasus dimana fungsi yang ditaksir diketahui bersifat periodik, maka nilai K1 dan 
K2 dihitung dengan rumus (4.114) dan (4.115) dan kemudian nilai f! —f!,, dihitung 
sebagai nilai rata-rata K1 dan K2, yaitu 

Ka 


Kw — HK, tK,) 


Selanjutnya nilai f," dan ft, dihitung sebagai berikut 








Data 3 (4.118) 
1 1 
1 fa H-f. K 
fu fu —ghora ora 3 av (4.119) 
YTA 2 8 Axa AXyra— 


Sekarang perhatikan bahwa persamaan (4.109) dapat ditulis secara simbolik menjadi 
bentuk berikut untuk 2 j« JTA —1 

ai kb hei 5 dj (4.120) 
dimana koefisien a, b, c dan d adalah sebagai berikut. 
Untuk 3S j JTA —2 


a Ang, 
bj - Aan, , An) 


Cc, 5 Ax 


d.-6 ra -f, Ran f Afi 
? Ax, Ax, 


j 1 


Untuk j-2 dan j-JTA-1 
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a,—0 


b, — (Ax, #Ax,) untuk kasus I 


b, - LA, #2Ax, untuk kasus II, III dan IV 








| untuk kasus I 


untuk kasus II 











#3K,, untuk kasus IV 


4 

AH 
#H 
Ha 


5 aan, untuk kasus III 


bra “2Axyrp ot Ata untuk kasus II, III dan IV 


























Cyra2 50 
Ka 6 tra — fora— | aa Ne 'ma| untuk kasus I 
Axa Axyra—3 
I f 
dya, 56 ora “ora | 4 aa “ma | untuk kasus II 
Ce Kxra— 1 Axyra- 2 
dya 56 fara “ora | | Kana aa oma tua 3K, untuk kasus III 
Axa 1 —2 
dia ha aa “oya | | Sama aa “ora 2 3K,, untuk kasus IV 
Axa 1 Kxyra 2 
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Sistem persamaan (4.120) adalah sebuah sistem persamaan aljabar linier simultan 
berbentuk tridiagonal dan dapat diselesaikan secara efisien dengan menggunakan 
algoritma Thomas. Prosedur metoda ini dapat dijelaskan sebagai berikut : 

Semua nilai koefisien a, b, c dan d untuk 2 £ j X JTA —1 telah dihitung dan diketahui. 


Prosedur eliminasi : 
Penetuan nilai awal 


B, —b, 
6, -—d, 
Secara beruntun hitung untuk j-3, 4, 5, ...., JTA-1 
k-aj/Bj 
Bj 5b —k.c,, 
6, —-d-k.8,, 
Prosedur Substitusi Balik 
Nilai awal : Era — ra Bora 
Secara beruntun dihitung untuk j-JTA-2, JTA-3,...., 4, 3, 2 
fp — (3, — #8, 
Setelah semua nilai fj" dihitung untuk 2 £ j X JTA —1 kemudian hitung nilai berikut 
f" 50 untuk kasus I 





f — Lg sah untuk kasus II 
1 
fa — 2 3— 5 2 -3 5, untuk kasus III 
Axi— Ax, 
f" Si 13— P aa, 3 — Kes untuk kasus IV 
Ar? “Ar, 


fara 50 untuk kasus I 


PL NN safara masi untuk kasus II 
2 JTA-1 


fa Te 





1 fr —f K 
bea Aa 3 Haa TA1 3 ——2 untuk kasus INI 
2 Ka X yA 
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For Hara -3 Kar untuk kasus IV 


X yra—1 x 





frra - Haa t3 
JTA- 


Kemudian nilai fj dihitung dengan rumus (4.96) untuk 1 X j £ JTA — 1 sebagai berikut. 


Kh “1 
1 4d u " 
fj ie maan ka Jan, 


Kasus untuk j-JTA nilai fj dihitung dengan rumus berikut 


SA MA | 
Fa - Tya Tarang 19 Ten al dna 


JTA-1 


Rumus interpolasi untuk selang j adalah 
£(x) — ((c1,.(x—x)-02).lx—x,) 1-03). Ja, (4.121) 
dimana koefisien C1, C2 dan C3 adalah 


Lt —f" 
C1 —— 
8 6 Ax, 

1 " 

C2, - 35 
C3, -fj 


Pengertian tentang variabel kompleks dapat diterapkan untuk menurunkan rumus 
interpolasi untuk memperbanyak jumlah titik pada permukaan aerofoil sembarang, 
dengan menggunakan fungsi interpolasi kubik lentur dalam variabel kompleks sebagai 
berikut : 


Data yang diberikan adalah koordinat (xAA,, YAA,) dari sejumlah JTA titik pada 


permukaan aerofoil ZA. Koordinat tersebut dapat diwakili oleh variabel kompleks ZAA 


dimana 
ZAA, — XAA, ti. YAA, (4.122) 


Sekarang kita tentukan bahwa aerofoil ZA dapat dipetakan menjadi lingkaran lex —1 di 
bidang-C oleh rumus z(C) sedemikian rupa sehingga peta dari ZAA j adalah GL, dimana 


UL, — expli.OL, ) - cos0L, &i.sin0L, (4.123) 
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